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PREFACE 


This  report  describes  the  logic  and  operation  of  the  simplified  spares 
optimization  model  called  DRAMA.  The  model  has  been  used  to  understand  better 
the  relationships  among  availability,  reliability,  and  spares  investment 
levels  for  new  weapon  systems. 
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EXECUTIVE  SUMMARY 


DRAMA  is  an  analytic  model  that  relates  the  availability  of  a  weapon 
system  to  the  cost  of  its  spares.  It  is  a  model  of  a  multi-item,  multi¬ 
location,  two-echelon  supply  system  in  which  each  item  may  be  stocked  at  a 
depot  and  two  or  more  similar  bases.  Item  demands  are  assumed  to  be  Poisson- 
distributed.  DRAMA  focuses  on  only  one  weapon  system  at  a  time,  so  users  who 
are  interested  in  modeling  multiple  weapon  systems  are  referred  to  the  LMI 
Aircraft  Availability  Model,  of  which  DRAMA  is  a  simplified  version. 

DRAMA  produces  a  set  of  points  comprising  an  availability-versus-cost 
curve.  Each  point  is  an  optimum,  i.e.,  it  represents  the  least-cost  mix  of 
spares  for  that  level  of  availability  or,  conversely,  the  greatest  avail¬ 
ability  for  that  level  of  spares  investment.  Availability  is  defined  as  the 
probability  that  an  end  item  (such  as  an  aircraft  or  tank)  is  not  waiting  for 
a  component  to  be  repaired  or  be  shipped  to  it;  thus,  availability  is  defined 
from  the  point  of  view  of  the  supply  system. 

The  potential  user  should  feel  encouraged  to  use  DRAMA  for  certain  ap¬ 
plications  for  which  it  may  not  appear,  at  first  glance,  to  be  entirely  ap¬ 
propriate.  For  example,  despite  its  assumption  of  a  two-echelon  supply  system 
in  which  all  bases  are  equal,  DRAMA  has  successfully  modeled  the  Army's  five- 
echelon  system,  and  has  represented  situations  where  bases  were  not 
equivalent. 

In  comparison  with  some  spares-optimization  models,  DRAMA  requires  rela¬ 
tively  little  input  data.  For  each  component,  DRAMA  needs  removal-rate, 
echelon-of-repair ,  cost,  and  application  data.  For  the  end  item,  it  needs 
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utilization  and  deployment  data  and  base  repair,  depot  repair,  order-and- 
shipment,  and  condemnation  lead  times. 

DRAMA's  primary  output  is  a  table  of  spares  costs  at  various  levels  of 
availability.  Optional  outputs,  including  quantities  of  items  and  costs,  are 
also  available. 
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I.  THE  ALGEBRAIC  MODEL 


THE  INVENTORY  SYSTEM 

DRAMA ^  is  a  model  of  a  multi-item,  multi-location,  two-echelon  inventory 
system  consisting  of  two  or  more  identical  using  locations  called  bases  and  a 
single,  central  resupply  location  called  a  depot.  Each  base  is  equipped  with 
end  items,  such  as  aircraft  or  tanks.  The  end  items  are  comprised  of  various 
components  critical  to  their  operation.  These  components  are  usually  called 
Line  Replaceable  Units  (LRUs)  and  are  to  be  distinguished  from  Shop  Replace¬ 
able  Units  (SRUs)  which  are  subcomponents  of  LRUs  and  which  are  not  directly 
removable  from  end  items.  DRAMA  models  only  LRUs,  not  SRUs. 

The  purpose  of  DRAMA  is  to  maximize  the  availability  of  end  items  subject 
to  a  constraint  on  the  total  cost  of  spare  components  or,  conversely,  to  find 
the  least  cost  mix  of  spare  components  that  will  result  in  a  specified  avail¬ 
ability.  Availability  is  defined  as  the  probability  that  an  end  item  selected 
at  random  will  not  be  waiting  for  a  component  to  be  repaired  or  be  shipped  to 
it. 

When  a  component  fails,  it  is  removed  from  the  end  item  and  sent  to 
base  repair.  If  a  serviceable  spare  component  is  available,  it  is  installed 
on  the  end  item.  (By  serviceable  we  mean  in  working  condition,  suitable  for 
use  in  the  end  item.)  The  base  repair  shop  decides  whether  the  repair  of  the 
failed  component  is  within  its  capability.  If  it  is,  it  repairs  the  component 
and  gives  it  to  base  supply  as  a  serviceable  spare. 

'''DRAMA  is  an  acronym  for  Diagnostic  Reliability,  Availability,  and 
Modularity  Analyzer.  The  name  comes  from  its  original  application,  which  was 
to  assess  the  impact  of  modularization  and  fault-diagnosis  reliability  on 
spares  costs. 
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If  the  repair  of  the  failed  component  is  beyond  the  repair  capability  of 
the  base,  it  is  declared  to  be  not-repairable-this-station  (NETS)  and  shipped 
to  depot  repair.  In  this  event,  a  replacement  item  is  immediately  ordered  by 
the  base  from  the  depot.  Thus,  the  inventory  system  operates  according  to  an 
(S-1,S)  inventory  policy,  i.e.,  whenever  the  number  of  components  on-hand  at 
the  base  plus  due-in  from  the  depot  minus  due-out  to  the  operational  organiza¬ 
tion  falls  below  the  authorized  stock  level ,  S,  another  component  is  ordered 
immediately. 

The  flow  of  components  in  the  system  is  described  in  Figure  1-1. 


figure  I  -  i 
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The  average  time  that  elapses  from  the  removal  of  a  failed  component  from 
an  end  item  until  it  is  returned  to  serviceable  condition  (given  that  it  is 
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not  NETS)  by  the  base  repair  shop  is  called  the  base  repair  time .  The  com¬ 
ponents  in  the  process  of  repair  at  the  base  are  said  to  be  in  the  base  repair 
pipeline .  We  define  a  pipeline  as  an  expected  value,  e.g.,  the  base  repair 
pipeline  of  a  component  is  the  expected  number  of  that  component  in  repair  at 
the  base. 

Similarly,  the  average  time  that  elapses  from  the  removal  of  a  failed 
component  from  an  end  item  (in  the  case  where  the  component  is  NETS)  until  it 
is  returned  to  serviceable  condition  by  the  depot  is  called  the  depot  repair 
time.  Note  that  the  depot  repair  time  includes  the  time  required  to  ship  the 
failed  component  from  the  base  to  the  depot.  That  shipment  is  referred  to  as 
retrograde  shipment.  Thus,  the  depot  repair  pipeline  includes  the  retrograde 
pipeline . 

The  average  time  that  elapses  from  the  ordering  of  a  serviceable  com¬ 
ponent  by  the  base  from  the  depot  until  the  component  is  received  by  the  base 
is  called  the  order-and-ship  time.  We  assume  that  the  determination  of  NETS 
status  and  the  placement  of  the  order  are  instantaneous. 

Some  proportion  of  the  NETS  components  sent  to  the  depot  for  repair  may 
also  be  beyond  the  repair  capability  of  the  depot,  or  the  cost  of  repair  may 
exceed  the  cost  of  a  replacement  item.  In  this  case,  the  depot  condemns  the 
component  and  orders  a  replacement  from  the  manufacturer.  The  average  time 
that  elapses  from  the  removal  of  such  a  component  from  an  end  item  until  it  is 
received  by  the  depot  from  the  manufacturer  is  called  the  condemnation  lead 
time  or  procurement  lead  time .  This  lead  time  is  combined  with  the  depot 
repair  time  for  our  purposes  so  that  whenever  we  refer  to  the  depot  repair 
pipeline,  we  include  the  condemnation  pipeline. 

Any  component  that  is  in  the  base  repair,  depot  repair  (including  retro¬ 
grade)  ,  or  order-and-ship  pipeline  is  said  to  be  in  resupply.  The  resupply 
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pipeline  (or  simply  pipeline)  is  the  sum  of  these  three  component  pipelines . 
The  average  time  a  component  is  in  the  resupply  pipeline  is  called  the  re¬ 
supply  time. 

The  inventory  of  spare  components  is  positioned  at  the  depot  as  depot 
stock  and  at  the  bases  as  base  stock .  Whenever  an  order  for  a  serviceable 
component  is  received  by  a  base  and  the  base  is  out  of  stock,  a  base  backorder 
occurs.  Similarly,  when  an  order  is  received  by  the  depot  and  the  depot  has 
no  stock,  a  depot  backorder  occurs.  The  time  that  elapses  from  the  receipt  of 
an  order  by  the  depot  until  a  serviceable  component  is  available  at  the  depot 
is  called  the  depot  delay  time .  If  one  or  more  serviceable  spares  are  in 
depot  stock  when  the  order  is  received,  the  depot  delay  time  is  zero. 

THE  MODEL 

The  fundamental  purpose  of  DRAMA  is  to  produce  a  set  of  points  describing 
the  relationship  between  end-item  availability  and  the  initial  investment  cost 
for  spare  components  for  depot  and  base  stocks.  Each  point  represents  the 
least-cost  mix  of  spares  for  its  corresponding  availability.  In  other  words, 
DRAMA  produces  a  spares  list  that  maximizes  end-item  availability  given  a 
budget  constraint;  conversely,  it  minimizes  spares  costs  for  any  specified 
level  of  end-item  availability. 

We  will  now  show  how  DRAMA  models  this  logistics  system.  We  assume  that 
demands  for  components  are  equivalent  to  removals  of  components  due  to 
failures  and  that:  (1)  demands  follow  a  Poisson  distribution  and  (2)  demands 
for  any  component  are  independent  of  demands  for  any  other  component.  Then, 
we  define  for  any  particular  component: 

T  =  average  resupply  time, 
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=  average  base  repair  time, 

=  average  depot  repair  time , 
t  =  average  order-and-ship  time, 

=  average  depot  delay  time, 
t^  =  average  condemnation  lead  time, 

<j)  =  fraction  of  NRTS  components  condemned  by  the  depot, 
r  =  base  repair  fraction  (i.e.,  1  -  NRTS  rate), 

A,  =  demand  rate, 

X  =  a  random  variable:  the  number  of  components  in  resupply,  and 
B(Sq,  S)  =  a  random  variable:  the  number  of  base-level  backorders  given 
a  depot  stock  level  of  Sq  and  a  base  stock  level  of  S. 

In  accordance  with  our  Poisson  assumption,  the  probability  density  func¬ 
tion  (p.d.f.)  of  X  is  given  by 


pOO  * 


e  AT  (XT)*,  x  =  0,  1,  2, 


x  i 


0  elsewhere 


where  T  is  a  weighted  combination  of  pipeline  times  given  by 

T  =  r  t,  +  (1  -  r)  (t  +  t  ). 
b  o  w' 

Despite  the  fact  that  T  is,  in  fact,  a  random  variable,  we  may  use  it  here 

2 

like  a  constant  due  to  a  theorem  and  proof  of  C.  Palm.  Treating  T  as  a  con¬ 
stant,  we  can  express  the  expected  number  of  components  in  resupply  as 
E(X)  =  A.T. 

The  expected  number  of  base-level  backorders  (EBO) ,  given  a  base  stock 

level  of  S,  and  a  depot  stock  level  of  Sq,  is  simply 

E[  B(S  S)  ]  =  I  p(x)  (x  -  S). 
x>S 

2 

Palm,  C., "Analysis  of  the  Erlang  Traffic  Formula  for  Busy-Signal 
Arrangements",  Ericsson  Technics,  No. 5,  1938,  pp. 39-58. 


1-5 


The  depot  delay  time  is  a  function  of  the  depot  stock  level  and  the  depot 
demand  rate.  The  greater  the  depot  stock  level,  the  less  the  average  waiting 
time  for  a  serviceable  spare.  Since  we  assume  a  steady-state  system,  the 
expected  depot  delay  time  is  equal  to  the  product  of  the  expected  number  of 
depot-level  backorders  and  the  mean  time  between  depot  demands,  i.e., 

E[D(So)] 

Cw  A(l-r) 

where  D,  a  random  variable,  denotes  the  number  of  depot-level  backorders  and 
A.  (1-r)  is  the  depot  demand  rate.  We  define  N  as  the  number  of  components  in 
the  depot  repair  pipeline  (including  the  condemnation  pipeline).  By  Palm's 
Theorem,  N  has  a  Poisson  distribution.  Its  expectation  is  equal  to 
A  =  A  (1  -  r)  [<J)tc  +  (1  -  <t>))td] 
and  its  p.d.f.  is  simply 


g(n) 


e~V\  n  =  0,  1,  2,  .  .  . 
n! 

0  elsewhere  . 


Finally,  E[D(S  )],  the  expected  number  of  depot-level  backorders  is  computed 
as 

EFD(S  )]  =  Z  g(n)  (n-S  ). 

°  n>S 

o 


We  will  now  extend  the  model  in  a  straightforward  manner  to  the  multi¬ 
item,  multi-base  system.  In  our  notation,  we  will  let  i  denote  the  ith  com¬ 
ponent  type  where  there  are  n  different  component  types  and  we  will  assume 
that  there  are  k  identical  bases.  Then,  for  component  i, 

=  a  random  variable:  the  number  of  components  of  type  i  in  resupply, 
T^  =  resupply  time, 
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t  .  =  depot  delay  time, 
w ,  1 


<t>^  =  fraction  of  NETS  components  condemned  by  the  depot, 

=  base  repair  fraction, 

=  demand  rate  at  each  base, 
p^(x)  =  the  p.d.f.  of  Xj, , 

S  .  =  depot  stock  level, 

=  stock  level  at  each  base, 

Bf(So  ^,S^)  =  a  random  variable:  base-level  backorders  at  each  base, 

D. (S  . )  =  a  random  variable:  depot-level  backorders, 

-L  U  j  1 

=  a  random  variable:  number  of  components  of  type  i  in  the  depot  repair 
pipeline,  and 

g.(n)  =  the  p.d.f.  of  N.. 

i  i 

The  variables  t^,  t^,  tQ,  and  t  are  defined  as  before;  they  are  assumed  to 
apply  to  all  components. 

Thus , 


n>S  .  x 

0,1 

and  E(Ni)  =  k  A  (1-r^  [<t>itc  +  (1  -  (J>i)td] . 
Denote  ECN^)  by  A..  Then, 


g,  (n)  =  • 


i 


n ! 


0  elsewhere  - 


The  depot  delay  time  for  component  i  is 


tw, i  kAi(l-ri) 
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The  p,d.f.  of  X.  becomes  simply 


pjU)  - 


-X.T. 

e  11  (X . T . ) x 
1  1 


X 


^  0  elsewhere 

where  T.  =  r.t,  +  (1  -  r.)  (t  +  t  . )  ; 

i  lb  i  v  o  w,i;  * 


so 


E[B.CS  .  ,S  )]  =  Z  p  (n)(n-Si). 
1  0,1  x>S . 

l 

THE  CONCEPT  OF  AVAILABILITY 


We  will  now  consider  the  relationship  between  base-level  backorders  and 
end-item  availability.  Suppose  component  i  has  a  quantity  per  end  item  of  q^, 
that  there  are  u  end  items  at  each  base,  and  that  there  are  k  bases  and  n 
different  component  types  in  the  system,  as  before. 

An  estimate  of  the  probability  that  component  i  is  missing  from  any  one 
of  its  positions  in  an  end  item  is 

E [B . (S  ,,S.)] 

1  0,11  . 


uq. 


The  estimated  probability  it  is  not  missing,  then,  is 

E [B . (S  . ,S.)] 

1  -  i  o,i  i 


uq, 


and  that  it  is  not  missing  in  any  of  its  positions  on  an  end  item  is 


E[Bi(S0>i.5i)] 


1  - 


^i 


qi 


Finally,  the  estimated  probability  that  no  component  is  missing  on  an  end  item 

q  • 


is 


n 

a  =  n 

i=l 


e[b.(s  .  ,S.)] 

u  0,1  3  J 


uqi 
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"A"  denotes  the  end-item  availability  rate.  It  is  this  function  that  DRAMA 
maximizes  subject  to  a  constraint  on  spares  investment. 

Maximizing  A  is  equivalent  to  maximizing  the  logarithm  of  A;  indeed,  that 
is  the  objective  function  operated  on  by  DRAMA.  This  objective  function  is 
separable ;  i.e.,  it  allows  us  to  compute  the  ratio  of  marginal  availability 
per  unit  cost  for  each  component  type  independent  of  all  other  component 
types.  The  actual  solution  technique  employs  marginal  analysis.  It  proceeds 
in  an  iterative  fashion.  At  each  iteration  it  computes  the  end-item  avail¬ 
ability  that  would  result  from  an  incremental  addition  of  spare  components  to 
the  system  for  each  component  type.  It  then  selects  the  component  type  that 
results  in  the  greatest  increase  in  end-item  availability  per  unit  cost. 
These  computations  are  based  on  the  allocation  of  spare  components  among  the 
depot  and  all  the  bases  that  yields  the  minimal  number  of  expected  base-level 
backorders  in  the  system.  The  solution  technique  is  explained  further  in 
Chapter  II. 

Three  fundamentally  important  assumptions  underlie  DRAMA'S  computations. 
The  first  is  that  component  demands  are  independent  of  one  another.  The 
second  is  that  the  process  that  generates  demands  (e.g.,  a  flying  hour  or 
sortie  program  for  aircraft  or  a  number  of  miles  of  operation  for  tanks)  is 
known  and  does  not  change  as  a  function  of  spares  levels.  The  third  assump¬ 
tion  is  that  the  system  is  in  a  steady  state,  i.e.,  the  demand  process  is 
stationary. 


1-9 


II.  THE  COMPUTER  IMPLEMENTATION 


STRUCTURE 

To  economize  on  core,  the  input  and  output  functions  of  DRAMA  have  been 
delegated  to  separate  programs.  These  programs  are  respectively  named  PROLOG 
and  EPILOG,  and  are  discussed  in  Chapter  III.  This  chapter  focuses 
exclusively  on  the  main  program,  DRAMA,  and  its  subroutines. 

DRAMA  consists  of  a  main  program  and  twelve  subroutines.  The  logical 
links  between  subroutines  are  shown  in  Figure  II - 1 .  The  functions  of  each 
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subroutine  are  discussed  briefly  below,  and  are  more  fully  described  in  the 
remaining  sections  of  this  chapter.  We  now  define  the  critical  variables  that 
the  subroutines  employ.  These  variables  are: 

BEBO:  The  expected  number  of  backorders  of  a  component  at  each  base. 

It  is  a  function  of  the  spares  inventories  at  the  depot  and  at 
each  base.  In  Chapter  I,  it  was  defined  as  B[S  ,S]. 

DEBO:  The  expected  number  of  backorders  of  a  component  at  the  depot. 

It  is  a  function  of  the  spares  inventory  at  the  depot,  and  it  was 
defined  in  Chapter  I  as  E[D(Sq)]. 

BMEAN:  The  mean  value  of  the  Poisson  distribution  used  in  the  computa¬ 
tion  of  BEBO.  BMEAN  is  equal  to  the  expected  number  of  back¬ 

orders  that  a  base  would  have  if  it  had  no  spares,  E[B(S  ,  0)]. 
Initially,  when  neither  the  bases  nor  the  depot  have  spares, 
BMEAN  is  also  equal  to  the  expected  number  of  parts  in  resupply 
to  each  base.  This  initial  value  of  BMEAN  is  a  distinct  variable 
in  DRAMA  called  BMEAN0.  In  Chapter  I's  notation,  BMEAN  is 
equivalent  to  XT. 

DMEAN:  The  mean  value  of  the  Poisson  distribution  used  in  the  computa¬ 
tion  of  DEBO.  It  is  equal  to  the  expected  number  of  parts  in  the 
depot  pipeline,  where  the  depot  pipeline  is  understood  to  include 
both  the  condemnation  pipeline  and  the  retrograde  pipeline,  in 
addition  to  the  depot  repair  pipeline.  DMEAN  is  equal  to  the 
expected  number  of  backorders,  DEBO,  when  the  depot  has  no 
spares.  In  Chapter  I's  notation,  DMEAN  is  equivalent  to  A. 
Unlike  BMEAN,  it  does  not  change  in  value  during  a  DRAMA  run. 
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The  subroutine  functions  can  be  summarized  as  follows: 

(1)  PRELIM:  This  is  a  preliminary  routine,  not  involved  at  all  in 
marginal  analysis.  At  the  user's  option,  it  will  "purchase"  a 
reasonable,  but  not  optimal,  quantity  of  each  component  for  the 
depot  and  for  each  base. 

(2)  AVAIL:  This  subroutine  computes  component  and  end-item  availabili¬ 
ties  . 

(3)  MARGIN:  This  routine  computes  a  marginal  benefit-to-cost  ratio  for 
each  component.  In  the  process  of  computing  these  ratios,  MARGIN 
determines  the  purchase  quantity  that  yields  the  maximum  ratio,  and 
it  computes  the  optimal  base/depot  distribution  of  the  quantity  of 
spares  that  would  be  available  after  the  purchase. 

(4)  COMPOP:  This  subroutine  ranks  the  components  according  to  their 
benefit-to-cost  ratios.  It  also  selects  the  component  with  the 
highest  ratio  when  a  purchase  has  to  be  made. 

(5)  UPDATE:  This  subroutine  recomputes  component  inventories,  system 
spares  costs,  and  component  and  end-item  availabilities  to  account 
for  a  spares  purchase. 

(6)  LIB:  Computes  expected  backorders  at  the  depot  (DEBO)  when  the 
depot  has  a  quantity  of  spares  equal  to  the  depot  pipeline  (DMEAN). 
Alternatively,  it  computes  expected  backorders  at  each  base  (BEBO) 
when  each  base  has  a  quantity  of  spares  equal  to  BMEAN. 

(7)  CON:  Computes  DEBO  given  a  quantity  of  depot  spares  equal  to  DMEAN 
minus  six  standard  deviations  (a  standard  deviation  equals  the 
square  root  of  the  mean).  Alternatively,  it  computes  BEBO  given  a 
quantity  of  spares  at  each  base  equal  to  BMEAN  minus  six  standard 
deviations . 
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(8)  DNEXT :  Computes  the  DEBOs  that  a  component  could  have  after  its 

next  spares  purchase. 

(9)  BNEXT :  Computes  the  BEBOs  that  a  component  could  have  after  its 

next  spares  purchase. 

(10)  SIZOP :  Determines  the  purchase-size  that  will  yield  the  greatest 
benefit  for  its  cost. 

(11)  FACTLN:  Computes  the  logarithm  of  the  probability  that  exactly  N 

units  will  be  in  the  depot  pipeline.  Alternatively,  it  computes  the 
logarithm  of  the  probability  that  exactly  N  units  would  be  back¬ 
ordered  at  each  base  if  it  had  no  spares. 

(12)  POISSY:  This  subroutine  computes  the  value  of  BEBO  that  would 

result  from  adding  K  spare  units  to  each  base. 

WALKTHROUGH 

DRAMA,  the  main  program,  initializes  all  of  the  variables  involved  in  the 
optimization  process.  For  each  component,  it  sets  the  counter  of  depot  spares 
(JSPARE)  and  the  counter  of  base  spares  (KSPARE)  to  zero.  In  addition,  it 
sets  BEBO  and  DEBO  equal  to  BMEAN0  and  DMEAN,  and  it  computes  values  for  all 
of  the  variables  that  are  later  involved  in  the  recomputation  of  BEBO  and 
DEBO.  These  variables  are  discussed  extensively  in  the  last  section  of  this 
chapter.  Finally,  DRAMA  sets  the  user-option  variables  (described  in  Chapter 
III)  to  their  user-chosen  values. 

One  of  the  options  at  the  user's  disposal  is  the  option  of  making  limited 
purchases  of  spares  prior  to  marginal  analysis.  This  is  done  by  calling 
PRELIM.  PRELIM  is  a  particularly  useful  option  for  the  user  who  is  interested 
in  observing  cost-availability  relationships  at  only  high  levels  of  avail¬ 
ability,  for  it  considerably  shortens  the  total  running  time  of  the  model 
without  impairing  its  accuracy.  However,  for  the  user  who  is  interested  in 
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observing  low  levels  of  availability,  PRELIM  is  less  desirable.  Because  it 
does  not  buy  "optimal"  quantities  of  each  component,  PRELIM  may  result  in 
excessive  costs  at  low  availability  levels  (at  higher  levels  of  availability, 
marginal  analysis  purchases  will  have  corrected  for  these  suboptimalities)  and 
because  running-time  problems  are  less  severe  at  these  levels,  the  advantages 
of  PRELIM  are  less  significant. 

If  the  user  wishes  to  exercise  his  option  of  making  a  preliminary  buy,  he 
must  assign  an  appropriate  value  to  the  user-option  variable  IPRE.  By  setting 
IPRE  equal  to  two  the  user  instructs  PRELIM  to  buy  liberally,  and  by  setting 
IPRE  to  one  he  instructs  the  subroutine  to  buy  conservatively.  Setting  IPRE 
to  zero  causes  DRAMA  to  bypass  PRELIM  entirely  and  to  proceed  directly  to 
marginal  analysis . 

When  instructed  to  buy  liberally,  PRELIM  buys  a  depot-pipeline's  worth  of 
spares  for  the  depot.  Then,  after  recomputing  the  expected  number  of  depot 
backorders  (DEBO)  and  the  Poisson-distribution  mean  at  each  base  (BMEAN)  ,  it 
buys  BMEAN  spares  for  each  base.  When  the  quantities  DMEAN  and  BMEAN  are  not 
integers,  PRELIM  buys  spares  quantities  equal  to  the  truncated  values  of  DMEAN 
and  BMEAN.  To  reinitialize  the  variables  associated  with  the  expected  back¬ 
order  computations,  PRELIM  must  call  the  subroutine  LIB.  This  process  is 
carried  out  for  each  component  in  the  end  item. 

When  instructed  to  buy  conservatively,  PRELIM  buys  a  quantity  of  spares 
for  the  depot  equal  to  DMEAN  minus  six  standard  deviations.  (Note:  the 
standard  deviation  equals  the  square  root  of  the  mean).  Then,  after  recomput¬ 
ing  the  values  of  DEBO  and  BMEAN,  it  buys  a  quantity  of  spares  for  each  base 
equal  to  BMEAN  minus  six  standard  deviations.  To  update  the  values  of  the 
variables  used  in  the  expected  backorder  computations,  it  calls  the  subroutine 
CON.  As  in  the  liberal  buy,  PRELIM  truncates  spares  purchase  quantities 
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whenever  BMEAN  and  DMEAN  are  non- integral ,  and  it  performs  its  computations 
independently  for  each  component. 

Whether  or  not  the  user  has  exercised  his  option  of  making  a  preliminary 
buy,  the  next  step  in  DRAMA  is  to  initialize  the  availability  of  every  com¬ 
ponent.  This  is  accomplished  by  calling  the  subroutine  AVAIL.  After  the 
availability  initialization,  DRAMA  begins  marginal  analysis. 

DRAMA's  first  step  in  marginal  analysis  is  to  compute  marginal  benefit- 
to-cost  ratios,  or  ranking  values,  for  each  component.  This  is  accomplished 
by  the  subroutine  MARGIN.  The  second  step  is  to  rank  the  components  according 
to  the  magnitude  of  their  ranking  values,  and  to  identify  the  component  with 
the  maximum  value.  This  is  the  task  of  the  subroutine  COMPOP.  The  third  step 
is  to  "buy”  spares  of  the  component  with  the  maximum  value.  This  is  ac¬ 
complished  by  the  subroutine  UPDATE.  Finally,  it  is  necessary  to  see  whether 
the  cumulative  spares  cost  or  the  system  availability  has  reached  a  con¬ 
straint.  If  either  of  the  constraints  has  been  reached,  marginal  analysis 
will  terminate.  If  not,  marginal  analysis  will  continue,  and  the  model  will 
return  to  step  one. 

In  the  second  and  subsequent  iterations  of  marginal  analysis,  it  is 
unnecessary  to  recompute’  every  component's  ranking  value.  This  is  because  the 
purchase  of  one  type  of  component  has  no  effect  whatsoever  on  the  ranking 
values  of  the  rest.1  Therefore,  when  MARGIN  is  recalled,  it  has  to  recompute 
only  one  ranking  value,  that  of  the  component  most  recently  purchased.  Simi¬ 
larly,  when  COMPOP  is  recalled,  it  merely  has  to  find  the  new  ranking  position 


The  benefit-to-cost  ratios  of  component  purchases  are  independent  of  one 
another  because  the  benefit  function  in  MARGIN  is  the  logarithm  of  end-item 
availability.  That  function,  unlike  the  availability  function,  is  separable 
by  component  type. 
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of  the  most  recently  purchased  component;  the  relative  positions  of  the  other 
components  remain  unchanged. 

When  marginal  analysis  ends,  DRAMA  outputs  a  set  of  up  to  27  cost-avail¬ 
ability  points  which  define  the  system's  optimal  cost-availability  curve. 
These  points  start  at  an  availability  level  of  five  per  cent  and  proceed  up  to 
an  availability  of  99  per  cent  as  follows:  5,  10,...,  90,  91,  92,...,  99 

per  cent.  If  the  availability  constraint  is  less  than  99  per  cent,  or  if  the 
cost  constraint  is  exceeded  before  99  per  cent  availability  is  achieved,  then 
only  a  partial  set  of  cost-availability  points  will  be  output. 

To  insure  that  the  quantity  of  spares  purchased  will  be  optimal,  and  to 
insure  that  each  component  has  an  optimal  distribution  of  spares  between  the 
depot  and  the  bases,  DRAMA  relies  on  the  optimizing  logic  of  MARGIN.  MARGIN 
calculates  a  component's  optimal  purchase  quantity  and  spares  distribution  in 
the  process  of  computing  its  benefit- to-cost  ratio;  then  it  returns  these  data 
to  DRAMA  so  that  the  component's  inventories  will  be  properly  updated  the  next 
time  it's  purchased.  Because  of  the  importance  of  MARGIN  to  the  overall 
process  of  optimizing  inventment  in  spares  across  components,  the  logic  of  the 
subroutine  is  discussed  extensively  below. 

THE  LOGIC  OF  MARGIN 

To  calculate  the  optimal  purchase  quantities  and  distributions  of  spares, 
MARGIN  engages  in  a  five-step  process.  The  steps  are: 

1)  Compute  the  value  of  DEBO  that  would  result  from  each 

purchase  alternative.  (A  purchase  alternative  is  a  particular 
purchase  quantity  distributed  in  a  particular  way.  There  are  a 
limited  number  of  purchase  alternatives  for  each  component.)  This 
is  done  by  DNEXT. 

2)  Compute  the  value  of  BMEAN  for  each  purchase  alternative.  (BMEAN  is 
a  function  of  DEBO.)  This  is  done  at  the  beginning  of  BNEXT. 
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3)  Compute  the  value  of  BEBO  for  each  purchase  alternative.  (BEBO  is  a 
function  of  BMEAN  and  the  base  spares  inventory) .  This  is  done  in 
the  latter  part  of  BNEXT,  and  it  may  involve  the  subroutines  LIB, 
CON,  FACTLN  and  POISSY. 

4)  Compute  the  benefit-to-cost  ratio  for  each  purchase  alternative 
(this  ratio  is  a  function  of  BEBO  and  the  purchase  quantity  involved 
in  the  purchase  alternative).  This  is  the  job  of  the  first  part  of 
COMPOP. 

5)  Find  the  maximum  ratio  and  identify  the  purchase  alternative  that 

yielded  it.  This  is  accomplished  in  the  latter  part  of  COMPOP. 

In  order  to  perform  these  computations  in  a  reasonable  amount  of  time, 
MARGIN  considers  only  a  limited  number  of  purchase  alternatives.  The  number 
of  alternatives  considered  can  never  exceed  the  number  of  bases  being  consid¬ 
ered  plus  one  (NOB+1).  Thus,  for  a  seven-base  scenario,  there  can  be  no  more 
than  eight  purchase  alternatives.  Those  alternatives  are  described  below  in 
Table  II-l,  with  j  denoting  the  current  spares  inventory  at  the  depot,  and  k 
denoting  the  current  spares  inventory  at  each  base. 

TABLE  II-l 

THE  PURCHASE  ALTERNATIVES  CONSIDERED 
FOR  A  SEVEN-BASE  SCENARIO 


PURCHASE 

ALTERNATIVE 

NUMBER 


QUANTITY  RESULTING  RESULTING 

PURCHASED  DEPOT  INVENTORY  BASE  INVENTORY 


1 

2 

3 

4 

5 

6 

7 

8 


1  j  -  6  k  +  1 

2  j  -  5  k  +  1 

3  j  -  4  k  +  1 

4  j  -  3  k  +  1 

5  j  -  2  k  +  1 

6  j  -  1  k  +  1 

7  j  k  +  1 

1  j  +  1  k 
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Note  that  MARGIN  needs  to  return  only  one  datum,  the  purchase-alternative 
number,  to  inform  DRAMA  of  the  optimal  purchase  quantity  and  the  optimal 
spares  distribution.  The  optimal  purchase  quantity  will  always  equal  the 
purchase-alternative  number  (NBEST)  unless  the  best  purchase-alternative  is 
alternative  eight,  in  which  case  the  purchase  quantity  will  equal  one.  The 
resulting  depot  inventory  will  always  equal  the  current  inventory  minus  the 
number  of  bases  plus  NBEST  (j  -  NOB  +  NBEST).  The  inventory  at  each  base  will 
equal  the  existing  inventory  (k)  plus  one,  unless  the  best  purchase  alter¬ 
native  is  alternative  eight.  Then  the  base  inventory  will  remain  the  same. 

Note  also  that  the  inventory  at  the  depot  can  decline  even  though  the 
total  number  of  spares  is  increasing.  This  is  because  higher  availabilities 
can  sometimes  be  achieved  simply  by  moving  spare  units  from  the  depot  to  the 
bases . 

The  set  of  permissible,  alternative  purchases  is  constrained  by  rules 
designed  to  minimize  an  item’s  expected  base-level  backorders.  Those  rules 
are : 

RULE  1 :  A  purchase  cannot  increase  the  inventory  at  any  site  by  more 
than  one  unit. 

At  any  site,  be  it  a  base  or  a  depot,  tjjLe  EBO  reduction  of 
additional  spares  is  strictly  decreasing.  The  second  spare 
reduces  EBO  less  than  the  first,  and  the  third  reduces  the  EBO 
less  than  the  second.  Therefore,  adding  one  spare  to  a  site 
will  always  be  more  cost-effective  than  adding  two  or  more. 

RULE  2 :  If  a  spare  is  added  to  the  depot,  none  can  be  added  to  the 
bases,  and  vice-versa. 

Consider  the  EBO-reduction- to-cost  ratios  of  the  following 
purchases:  the  purchase  of  a  single  depot  spare,  which  reduces 


Decreasing  marginal  EBO  is  always  present,  but  it  will  not  necessarily 
be  evident  until  a  site’s  inventory  of  an  LRU  exceeds  the  site’s  Poisson 
distribution  mean  minus  six  standard  deviations. 
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EBO  across  all  bases  by  the  quantity  D;  the  purchase  of  N  base 
spares  (no  more  than  one  per  base) ,  where  each  spare  reduces 
the  EBO  across  all  bases  by  the  quantity  B;  the  purchase  of  one 
depot  spare  together  with  N  base  spares,  where  the  depot  spare 
reduces  EBO  by  the  quantity  D,  and  where  each  base  spare  re¬ 
duces  EBO  by  the  quantity  B’ .  Because  the  EBO-reduction  of  a 
spare  decreases  as  BMEAN  decreases,  and  because  BMEAN  decreases 
as  the  depot  spares  inventory  increases,  B  is  greater  than  B' . 

To  compute  the  EBO-reduction-to-cost  ratios,  we  assume  for 
simplicity  that  each  spare  costs  $1.00.  Therefore,  the  EBO- 
reduction-to-cost  ratios  for  each  purchase  are: 

D  for  the  purchase  of  a  single  depot  spare, 

B  for  the  purchase  of  N  base  spares, 

and 

D+NB 1  for  the  purchase  of  a  depot  spare  together  with 
N+l  N  base  spares. 

Now,  if  the  depot  spare  is  at  least  as  cost-effective  as  the  N 
base  spares,  it  is  more  cost-effective  than  the  combined 
purchase : 


D  ^  B  >  B' 

ND  >  NB' 

(N+l)  D  >  D  +  NB' 

D  >  D  +  NB' 

N  +  1 

Similarly,  if  the  N  base  spares  are  at  least  as  cost-effective 
as  a  single  depot  spare,  then  they  are  more  cost-effective  than 
the  combined  purchase. 

B  ^  D 
NB  >  NB' 

(N+l)  B  >  D+NB' 

B  >  D+NB' 

(N+l) 

Therefore,  it  will  always  be  more  cost-effective  to  buy  either 
one  depot  spare  or  N  base  spares  than  to  buy  a  depot  spare 
together  with  N  base  spares.  As  a  consequence  of  this  rule, 
the  maximum  purchase  quantity  cannot  exceed  the  number  of 
bases . 

RULE  3:  All  bases  must  have  identical  spares  inventories. 

This  is  an  integral  part  of  the  model's  assumption  that  all 
bases  are  identical  in  every  respect. 
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RULE  4: 


If  the  purchase  quantity  is  greater  than  one,  but  less  than 

NOB,  spares  must  be  redistributed  from  the  depot  to  the  bases. 

Suppose  that  more  than  one  spare  is  purchased.  Rule  one  pro¬ 
hibits  us  from  putting  the  entire  purchase  at  the  depot,  and 
rule  three  prohibits  us  from  simply  putting  the  entire  purchase 
at  the  bases,  unless  the  purchase  quantity  is  NOB.  To  satisfy 
rule  three,  we  have  to  put  more  spares  at  the  bases  than  we 
purchased.  Therefore,  we  have  to  transfer  a  quantity  of  spares 
from  the  depot  to  the  bases.  That  quantity  is  equal  to  the 
number  of  bases  (NOB)  minus  the  number  of  spares  that  we  pur¬ 
chased  . 

RULE  5 :  In  no  case  will  more  than  NOB-1  spares  be  redistributed . 

We  have  not  shown  that  it  will  never  be  optimal  to  redistribute 
more  depot  spares;  however,  we  have  not  observed  any  cases  in 
which  greater  redistribution  will  result  in  lower  levels  of 
EBO. 

Applying  these  rules,  it  is  clear  that  there  can  only  be  NOB+1  per¬ 
missible  purchase  alternatives.  For  a  purchase  quantity  of  one,  there  are 
only  two  distributional  alternatives--the  one  spare  can  be  placed  at  the 
depot,  or  the  one  spare  can  be  placed  at  the  base,  with  NOB-1  spares  being 
redistributed.  And  for  purchase  quantities  of  more  than  one,  there  is  just 
one  distributional  alternative--all  of  the  purchased  spares  must  go  to  the 
bases,  and  NOB  minus  the  purchase  quantity  must  be  redistributed. 

COMPUTING  EXPECTED  BACKORDERS 


Most  of  the  time,  expected  backorders  at  the  base  and  the  depot  are  com¬ 
puted  with  an  iterative  algorithm.  We  describe  that  algorithm  below  as  it  is 
used  in  computing  expected  backorders  at  the  depot.  The  same  algorithm  is 
used  in  computing  the  expected  backorders  at  the  bases,  the  only  difference 
being  the  names  of  the  variables. 

The  algorithm  that  computes  DEBO  performs  the  following  steps  with  each 
iteration: 

STEP  1:  Compute  the  probability  that  the  population  in  the  depot  pipe¬ 

line  is  equal  to  the  inventory  on  hand.  The  logarithm  of  that 
probability  is  called  DPROB. 
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STEP  2:  Compute  the  reduction  in  backorders  that  will  result  from  the 

addition  of  a  spare.  This  quantity  is  a  function  of  DPROB,  and 
is  called  DREBO. 

STEP  3:  Compute  the  level  of  depot  backorders  that  will  result  from  the 

addition  of  a  spare.  This  variable  is  called  DEBO.  The  new 
value  of  DEBO  will  equal  the  existing  value  of  DEBO  minus 
DREBO. 

In  the  discussion  below,  we  will  denote  the  existing  value  of  DEBO  as 
DEBOjj ,  and  we  will  denote  the  new  value  as  DEBO^+^.  The  same  notation  will 
also  be  used  for  DREBO  and  DPROB. 

All  of  these  variables  can  be  expressed  in  forms  that  lend  themselves  to 
iterative  computations.  Consider  DPROB,  the  logarithm  of  a  Poisson  prob¬ 
ability.  We  know  that  the  probability  of  there  being  N  items  in  the  pipeline 
is 


-  DMEAN 
e 


DMEANN 


N! 


Thus  DPROB. T  is 


and  DPROB^+j  is  equal  to 


or,  equivalently, 


-DMEAN  +  N  In (DMEAN)  -  ln(N!) 

-DMEAN  +  (N  +  1)  ln(DMEAN)  -  ln[(N+l)!] 
DPROBN  +  In (DMEAN)  -  ln(N+l) 


DREBO  can  also  be  computed  iteratively.  To  see  how  it  can  be  computed 
iteratively,  consider  a  situation  in  which  there  are  no  depot  spares.  A  spare 
at  this  point  will  reduce  the  number  of  depot  backorders  to  zero  when  it 
otherwise  would  have  been  one,  to  one  when  it  otherwise  would  have  been  two, 


11-12 


to  two  when  it  otherwise  would  have  been  three,  etc.  In  short,  it  reduces 
DEBO  from 


03 

2  P.  (i) 

i=o 

to 

03 

I  P,  (i  -  1) 
i=l 

for  a  net  reduction  of 


CO 

DREBCL  *  E  p‘,  =  1— p 

,  1  •  T  1  0 

1=1 

where  P.  represents  the  probability  that  i  units  are  in 
the  deport  pipeline. 

Adding  a  second  spare  will  further  reduce  DEBO.  When  the  pipeline  contains 
two  units,  there  will  now  be  no  backorders;  when  it  contains  three  units, 
there  will  now  be  only  one  backorder,  etc.  The  net  reduction  in  DEBO  result¬ 
ing  from  the  second  spare  can  thus  be  expressed  as 

03 

DREB0o  =  I  p.  =  1  -  p  -  p  =  DREBO,  -  Pl 
2  ._0  ri  ro  rl  1  F1 

i-2 

Similarly,  the  net  reduction  in  DEBO  resulting  from  the  N  +  1st  spare  can 
be  expressed  as 

DREBOn+1  1  ■  P0  ■  Pi  '  P2  "  p3  •  •  •  PN  =  drebon  '  Pjf 


Since  is  simply  the  antilogarithm  of  DPROBN,  DREBO^+1  can  be  easily  com- 
puted . 

Using  the  DRAMA  variable  names,  the  iterative  computation  of  DEBO  can  be 
summarized  as  follows: 

DPROBn  =  DPROBN  +  ln(DMEAN)  -  In  (N) 
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DREBO„  =  DREBOXT  -  eDPR0BN 

N+l  N 

DEBO..  =  DEBO„  -  DREBOXTj, 

N+l  N  N+l 

The  appropriate  initial  values  of  these  variables  are  as  follows: 


DPROB 

=  -  DMEAN 

0 

DREBO 

=  1.0 

0 

DEBO 

=  DMEAN 

For  the  bases  a  similar  computation  of  EBO  is  carried  out  in  subroutine 
POISSY.  POISSY's  computations  can  be  summarized  as  follows  for  any  specified 
number  of  depot  spares: 

BPROBn  =  BPR0BN-1  +  ln(BMEAN)  -  In  (N) 

BREBOn+1  =  BEEBOj,  -eBPR0BN 

BEB0N+1  *  BEB0H  -  breb°!J+1 
The  appropriate  initial  values  of  these  variables  are 


BPROB 

0 

=  -  BMEAN 

BREBO 

=  1.0 

0 

BEBO 

=  BMEAN 

To  save  on  running  time,  the  model  sometimes  uses  shortcuts  instead  of 
this  iterative  computation.  One  shortcut,  the  subroutine  CON,  computes  BPROB , 
BREBO,  and  BEBO  for  an  inventory  of  spares  equal  to  BMEAN  minus  six  standard 
deviations;  the  other  shortcut,  the  subroutine  LIB,  computes  BPROB,  BREBO,  and 
BEBO  for  an  inventory  of  spares  equal  to  BMEAN.  These  routines  are  regularly 
called  by  the  subroutine  BNEXT.  They  may  also  be  called  by  PRELIM,  in  which 
case  they  will  compute  the  depot  variables,  DPROB,  DREBO,  and  DEBO,  as  well  as 
the  base  variables. 

The  logic  in  CON  is  based  on  the  fact  that  Poisson  probabilities  are 
infinitesimal  (less  than  10  when  N  is  less  than  BMEAN  minus  six  standard 
deviations  (BMEAN  minus  six  standard  deviations  is  called  ISIG) .  Because  of 
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this  fact,  the  first  ISIG  spares  each  provide  a  backorder  reduction  of  ap¬ 
proximately  one  unit,  and  the  value  of  BEBO,  after  their  purchase,  is  simply 
BMEAN  -  ISIG  . 

When  the  spares  inventory  is  greater  than  ISIG,  but  less  than  BMEAN,  CON 
has  to  compute  BREB0ISI(,  and  BPROBjgjg.  The  value  of  BREBO^g^  is  one,  be¬ 
cause  the  Poisson  probabilities  have  all  been  rounded  to  zero.  The  value  of 
BPROBjgjg  is  given  by  the  equation 

BPROBisig  =  -BMEAN  +  ISIG  ln(BMEAN)  -  STIRL 
where  STIRL  is  equal  to  the  ln(ISIG!)  and  is  computed  by  the  subroutine  FACTLN 
according  to  Stirling's  formula.  CON  then  transfers  these  values  to  POISSY, 
which  uses  them  to  iteratively  compute  BEBO^.  By  calling  CON  prior  to  POISSY, 
far  fewer  iterations  are  required  to  compute  BEBO^.  This  saves  a  great  deal 
of  time,  particularly  when  the  mean  of  the  base  Poisson  distribution  is  large. 

The  subroutine  LIB  is  based  on  the  equivalence  between  the  Poisson 
distribution  function  and  the  incomplete  gamma  function.  Because  of  this 
equivalence,  and  because  BREBO^+^  can  be  defined  as  the  complement  of  the 
Poisson  distribution  function,  BREBO^+1  is  equal  to  the  complement  of  the 
incomplete  gamma  function  evaluated  at  N.  According  to  a  derivation  published 

4 

by  Knuth,  the  complement  of  the  incomplete  gamma  function  is  approximately 


y-2/3  ' 

-1/2 

1 

'23  Y  Y3  " 

.V  211 

X 

u  270  12  6  . 

where  x  is  equal  to  BMEAN,  and 
where  y  is  equal  to  BMEAN  minus  the 
inventory  level. 


Knuth,  Donald  E.,  Fundamental  Algorithms,  Volume  I,  p.116,  Addison- 
Wesley  Inc.,  Reading,  Mass.,  1973. 
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Although  this  equation  can  approximate  BREBO  for  various  inventory 
levels,  the  approximations  are  most  accurate  when  the  inventory  level  is  equal 
to  the  truncated  value  of  BMEAN.  This  is  the  inventory  level  that  is  used  in 
LIB.  In  the  model,  this  inventory  level  is  called  IMEAN . 

LIB  computes  BPR0BIMEAN  just  as  CON  computes  BPROBISIG,  by  calling  FACTLN 
and  applying  the  appropriate  equation.  It  then  computes  BEBO^^  with  the 
equation  below: 

BPROB 

BEB0 IMEAN  =  (BMEAN- IMEAN)  (1-BREBO^^)  +  BMEAN  •  e  IMEAN  . 

When  the  stockage  level  at  each  base  exceeds  IMEAN,  DRAMA  calls  POISSY  to 
compute  the  backorder  reduction  of  the  remaining  stock. 
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III.  USING  THE  MODEL 


INPUTS 


The  DRAMA  main  program  requires  system-level  data  such  as  the  number  of 
bases  and  the  number  of  end  items  at  each  base,  and  some  component-level  data, 
such  as  cost  and  the  means  of  the  depot  and  base  Poisson  distributions,  DMEAN 
and  BMEAN0 .  Since  these  last  two  data  elements  are  not  usually  available  in 
weapon  system  data  bases ,  it  may  be  necessary  to  compute  them  from  elements 
that  are  available.  These  computations  are  carried  out  by  the  preprocessing 
program  called  PROLOG.  In  the  sections  below,  we  discuss  PROLOG'S  outputs  and 
the  manner  in  which  PROLOG  computes  DMEAN  and  BMEAN0. 


PROLOG 


The  standard  PROLOG  output  consists  of  scenario  data  describing  the 
number  of  bases  and  end  items  per  base,  followed  by  a  data  line  for  each 
component  type.  A  brief  description  of  the  necessary  data,  along  with  the 
corresponding  FORTRAN  variable  name,  is  given  below.  In  addition,  a  sample 
PROLOG  output  appears  in  Table  I I I — 1 . 

Scenario  Data:  NOB,  SPB,  NCOMP 


NOB 


Number  of  bases  supported  by  the  depot. 


SPB 


Number  of  end  items  per  base. 


NCOMP  -  Number  of  components . 

Component  Data:  I,  COST (I) ,  IPS (I) ,  DMEAN(I) ,  BMEAN0(I) 


I 


Counter  used  to  indicate  type-I  component. 


ICOST(I)-  Type-I  unit  cost. 


IPS(I) 


Installations  per  system;  number  of  type-I  components 


installed  on  each  end  item. 
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DMEAN(I) 


The  expected  number  of  backorders  at  the  depot  given 
zero  spares  of  type  I. 

BMEAN0(I)  -  The  expected  backorders  at  a  base  given  zero  spares 
of  type  I . 


TABLE  III-l.  SAMPLE  PROLOG  OUTPUT  (DRAMA  INPUT) 


17 

121. .0  70 

) Scenario  Data : 

)  NOB,  SPB,  NCOMP 

1 

1607.70 

1 

2.18 

8.12 

2 

1435 

6 

9.59 

105.53 

3 

227 

4 

326.97 

4591.79 

4 

1514 

1 

.96 

3.74 

Component  Data : 

5 

1514 

1 

.96 

3.74 

I,  ICOST(I) ,  IPS(I), 
BMEAN0(I),  DMEAN(I) 

6 

2230 

4 

2.00 

7.78 

7 

2230 

4 

2.00 

7.78 

8 

318 

2 

1.77 

6.87 

Scenario  Data :  SYSRAT ,  NSUBS 


SYSRAT  -  System  utilization  rate  for  the  particular  scenario  being 
modeled . 

NSUBS  -  The  number  of  subsystems  in  the  data  base. 

Component  Data:  RATE,  DRCT,  BRCT,  PLT,  OST 

RATE (I)  -  Component  removal  rate  per  application. 

ICND(I)  -  The  percentages  of  removed  units  that  get  condemned. 
IRTS(I,4)  -  Repairable-this-station  rates. 

DRCT  -  Depot  repair  time. 

BRCT  -  Base  repair  time. 
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DSRCT 

In  four-echelon  systems,  the  repair  time  at  the 

just  above  the  base  echelon. 

echelon 

GSRCT 

In  four-echelon  systems,  the  repair  time  at  the 

just  below  the  depot  echelon. 

echelon 

PLT 

Procurement  lead  time;  the  time  required  to 

a  condemned  unit. 

replace 

DOST 

Order-and-ship  time  from  the  depot  to  the  bases. 

GSOST 

In  four-echelon  systems,  the  order-and-ship  time 

echelon  just  below  the  depot  to  the  bases. 

from  the 

ISS(I)  - 

The  subystem  in  which  unit  I  is  installed. 

K(ISS)  - 

The  ratio  of  expected  subsystem  removal  rate  to  the  input 

subsystem  removal  rate. 

An  input  data  base 

for  PROLOG  is  illustrated  below. 

TABLE  III-2 .  SAMPLE  PROLOG  INPUT 


17 

1.0 

1.0 

121 

2.0 

70 

23.2 

i 

Scenario  data: 

NOB,  SPB,  NSUBS,NCOMP, SYSRAT 
Subsystem  data:  K(ISS) 

1607.00 

1 

10 

20 

30  40 

0  1  7.5 

Component  data: 

1435.00 

6 

5 

10 

20  40 

25  1  22.3 

ICOST(I) ,  IPS(I),  IRTS(I , 1) , 

• 

• 

• 

* 

* 

♦  •  ♦ 

IRTS (I , 2) ,  IRTS(I,3) , 

• 

• 

• 

• 

• 

♦  ♦  ♦ 

♦  ♦  ♦ 

IRTS(I,4) ,  CND(I),  ISS(I),  RATE ( I ) 

PROLOG 

computes 

the  rate 

of  removal  of 

a  component  at  an  individual  base 

as  follows: 

BRATE(I)  =  {K(ISS(I)-RATE(I)  }  {SYSRAT}  (IPS(I)-SPB} 

This  quantity  is  used  to  compute  the  quantities  BMEAN0  and  DMEAN. 


1 1 1 -3 


A  component's  depot  pipeline  (DMEAN(I))  consists  of  its  depot  repair 
pipeline,  its  condemnation  pipeline  and,  in  cases  where  four  echelons  are 
modeled,  its  "general  support"  pipeline.  Since  a  pipeline  has  a  population 
equal  to  the  rate  at  which  parts  enter  it  times  the  length  of  time  they  remain 
within  it,  these  three  pipelines  are  respectively  equal  to 
DPIPE(I)  =  IRTS (I ,4)  •  NOB  •  BRATE(I)  •  DRCT, 

CPIPE(I)  =  ICND(I)  •  NOB  •  BRATE(I)  •  PIT, 
and  GPIPE(I)  =  IRTS(I,3)  •  NOB  •  BRATE(I)  •  GRCT. 

The  depot  pipeline  equals  the  sum  of  these  three  pipelines,  namely: 

DMEAN(I)  =  DPIPE(I)  +  CPIPE(I)  +  GPIPE(I). 

A  component's  resupply  pipeline  at  an  individual  base,  BMEAN0 ,  consists 
of  a  base's  share  of  the  depot  pipeline  plus  the  base  repair,  direct-support 
repair,  and  order-and-shipping  pipelines.  These  pipelines  are  respectively 
BRPIPE  =  IRTS (1,1)  ♦  BRATE(I)  •  BRCT 
DRPIPE  =  IRTS (I, 2)  •  BRATE(I)  •  DSRCT 
GOPIPE  =  IRTS (I, 3)  •  BRATE(I)  •  GOST 
and  DOPIPE  =  [IRTS(I,4)  +  ICND(I)]*  BRATE(I)  •  DOST. 

The  resupply  pipeline  at  an  individual  base  is  given  by 

BMEAN0(I)  =  DPIPE('1^  +  BRPIPE  (I)  +  DRPIPE(I)  +  GOPIPE(IP)  +  DOPIPE(I). 

NOB 

STANDARD  DRAMA  OUTPUT 

Every  DRAMA  run  will  output  a  27-line  cost-availability  table.  A  sample 
output  appears  in  Table  III-3. 

OPTIONAL  DRAMA  OUTPUT 

DRAMA  will  also  output  special  files  at  the  user's  request.  (The  user 
must  attach  these  extra  files  to  the  DRAMA  program  in  order  to  output  data  to 
them.  The  appendix  explains  how  this  can  be  done  on  the  Honeywell  635)  . 
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TABLE  III-3.  SAMPLE  AVAILABILITY-COST  RELATIONSHIP 


Line 

No. 

Availability 

Total 

Spares 

Cost 

Line 

No. 

Availability 

Total 

Spares 

Cost 

1 

5% 

386061696 

15 

75% 

584437472 

2 

10% 

387825984 

16 

80% 

614785536 

3 

15% 

389122432 

17 

85% 

655098880 

4 

20% 

393912320 

18 

90% 

707731904 

5 

25% 

396275840 

19 

91% 

719241698 

6 

30% 

399317952 

20 

92% 

731577216 

7 

35% 

402902944 

21 

93% 

744964448 

8 

40% 

408072640 

22 

94% 

759649216 

9 

45% 

417048192 

23 

95% 

776184256 

10 

50% 

441347328 

24 

96% 

795014144 

11 

55% 

471371392 

25 

97% 

817839264 

12 

60% 

500581728 

26 

98% 

847312960 

13 

65% 

528524608 

27 

99% 

892400672 

14 

70% 

556713632 

DRAMA  can  output  four  special  files.  One  is  the  Base  Spares  File,  which 
contains  the  quantity  of  each  LRU  in  the  base  inventory  for  each  of  the 
27  availability  levels.  The  second  is  the  Depot  Spares  File ,  which  contains 
the  quantity  of  each  LRU  in  the  depot  inventory  for  each  of  the  27  availa¬ 
bility  levels.  The  third  is  the  BEBO  File,  which  contains  the  expected  number 
of  missing  units  for  each  LRU  for  each  of  the  27  availability  levels,  and  the 
last  is  the  DEBO  File,  which  contains  the  expected  number  of  outstanding 
orders  at  the  depot  for  each  LRU. 

To  obtain  these  files,  the  user  must  specify  an  output  option  in  the 
OPTION  section  of  the  main  program.  If  he  specifies  no  option,  he  gets  only 
the  standard  output.  If  he  sets  the  variable  IOUT  equal  to  one,  he  will  get 
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only  the  Base  Spares  File  in  addition  to  the  standard  output.  If  he  sets  IOUT 
equal  to  two,  he'll  get  both  the  Base  Spares  File  and  the  Depot  Spares  File. 
If  he  sets  IOUT  equal  to  three,  he'll  get  all  of  the  extra  files. 

EPILOG 

EPILOG  processes  DRAMA'S  output  files  and  generates  detailed  tables  of 
LRU  costs  and  backorders.  To  run  EPILOG,  the  user  must  input  the  DRAMA  output 
files  and  the  DRAMA  input  file  (The  method  for  inputting  these  files  is 
described  in  the  next  chapter.)  In  addition,  he  must  specify  his  EPILOG 
output  options. 

The  variable  IOUT  appears  again  in  EPILOG  and,  again,  its  purpose  is  to 
dictate  what  will  be  output.  By  setting  it  equal  to  one,  the  user  will  obtain 
a  listing  of  the  Base  Spares  File;  by  setting  it  equal  to  two,  the  user  will 
obtain  a  listing  of  the  Depot  Spares  File  as  well;  by  setting  it  equal  to 
three,  he'll  obtain  a  comprehensive  output  table  of  quantities  of  component 
spares  at  the  base  and  the  depot,  the  cost  of  the  component  spares  at  the  base 
and  the  depot,  the  unsatisfied  component  demand  (expected  back  orders  or  EBOs) 
at  the  base  and  the  depot,  and  the  component  availability  at  the  base.  A 
sample  EPILOG  output  for  IOUT  =  3  is  shown  in  Table  1 1 1 — 5 .  Notice  that  this 
output  is  for  a  special  case  where  there  were  initially  no  depot  backorders  of 
any  component. 

The  first  column  of  this  table  is  the  LRU  reference  number;  the  second  is 
the  quantity  of  spares  at  the  depot;  the  third  is  the  quantity  of  spares  at 
each  base;  the  fourth  is  the  total  spares  in  the  system;  the  fifth  is  the  cost 
of  the  depot  spares;  the  sixth  is  the  cost  of  the  base  spares;  the  seventh  is 
the  total  cost;  the  eighth  is  the  LRU's  DEBO;  the  ninth  is  the  LRU's  BEBO  and 
the  tenth  is  the  expected  LRU  availability.  In  addition  to  the  table,  EPILOG 
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will  output  system-level  summaries  of  base  and  depot  inventories  and  costs. 
The  summary  for  the  run  that  produced  Table  III-5  appears  in  Table  III-4, 
below. 


TABLE  III-4.  SUMMARY  STATISTICS  AT  THE  SYSTEM  LEVEL 


Total  Spares  Purchased 

894 

Total  Spares  at  the  Depot 

0 

Total  Spares  at  Each  Base 

298 

Total  Cost  of  All  Initial  Spares 

5773752 

Total  Cost  of  Spares  Bought  for 

the  Depot 

0 

Total  Cost  of  Spares  Bought  for 

the  Bases 

5773752 

Base  Expected  Back  Orders 

0.01489 

Depot  Expected  Back  Orders 

0 

Full  Systems  Availability 

99.5% 

Normally,  EPILOG  will  output  a  table  for  each  cost-availability  pair  in 
DRAMA's  standard  output.  Should  the  user  wish  to  obtain  an  output  table  for 
just  one  availability  level  he  can  do  so  by  setting  NAVAIL  equal  to  the  line 
number  (from  DRAMA's  standard  output)  of  the  availability  that  he  desires  to 
see.  The  line-numbers  of  each  availability  are  shown  in  Table  III-3.  Failure 
to  assign  a  value  to  NAVAIL  will  result  in  the  output  of  tables  for  all  of  the 
27  availability  levels. 
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TABLE  III-4 .  SAMPLE  PROLOGUE  OUTPUT  FOR  "TOUT11  EQUALS  3 


*  VVfM 

*  *  ¥  *  * 

*  *  *  *  * 

**  **  *******  sum*  Am  y  oi 

MODULE  PJkCHASES 

UHHHf  »*»***♦**  »****♦*♦*♦ 

*  ♦ 

»  **  ♦  *  * 

»  *  * 

♦  J  ^ 

A  215 

MODULE  SYSTEM 

rfITH  A  01A  GNOSTIC 

ACCUkACY  OF 

1QQ  *2  *********** 

*  * 

M  OU 

DEPO 

6m  Sc 

1 J  T 

£  OLPOr  $ 

%  OASES  % 

$*»  TOTAL  S$ 

£U  O- DEPOT 

EtiO-BASES 

l. 

J 

2 

b 

0  . 

11 11 2. 

1  111  2. 

0  . 

0. 00007366 

2 

u 

1 

3 

G  . 

38  95  8. 

3895  8. 

0. 

0. 00001129 

3 

li 

1 

3 

0. 

4  50  0. 

4500. 

0  . 

0.00000040 

4 

0 

•> 

b 

u . 

26  2b 8 . 

2  6  26  8  . 

0. 

0.00001071 

5 

0 

1 

3 

r . 

4  1  5  68  3  . 

41568  3. 

o . 

0  .  0  00  2i  364 

6 

c 

2 

6 

Q  . 

3082a. 

3  0  82  8. 

o  . 

0. 00000421 

7 

Li 

1 

3 

c . 

154  47  3. 

15447  3. 

0  . 

0.00161703 

a 

u 

2 

6 

• 

5283  6. 

52836. 

0  . 

0.00002683 

9 

0 

1 

3 
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GLOSSARY  OF  DRAMA  VARIABLES 

AM 

la  LIB  or  CON,  the  mean  of  Poisson  distribution  of  backorders 
at  either  the  depot  or  the  bases. 

AVCOMP 

The  availability  of  a  particular  component  in  EPILOG. 

AVLOG(I) 

The  logarithm  of  a  component’s  availability. 

AVTEST 

The  next  availability  milestone  at  which  availability  and 
cost  data  will  be  output.  Normally,  it  is  initialized  at 
five  per  cent. 

AVTOT 

A  particular  availability  level  in  EPILOG. 

BBOLIB(N) 

In  MARGIN,  the  value  that  BEBO  will  assume  if  purchase  alter¬ 
native  N  is  selected. 

BBOTOT 

BEBO  summed  across  all  components  at  a  particular  avail¬ 
ability  level.  Used  only  in  EPILOG. 

BOOST 

The  total  cost  of  base  spares  at  a  particular  availability 
level.  Used  only  in  EPILOG. 

BEBO(I) 

The  expected  number  of  backorders  of  component  (I)  at  each 
base. 

BMIOG 

In  BNEXT,  the  logarithm  of  the  variable  BMEANN,  which  denotes 
the  mean  value  of  the  base  Poisson  distribution  for  purchase 
alternative  N. 

BMEAN0(I) 

The  expected  number  of  units  of  component  (I)  in  resupply  to 
each  base.  Initially,  when  a  component  has  no  spares,  it  is 
also  equal  to  BEBO(I).  It  is  called  a  ’’mean”  because  it  is 
the  mean  value  of  the  Poisson  distribution  of  backorders  at 
each  base  when  the  depot  has  zero  spares. 

BMEANN 

In  MARGIN,  the  mean  of  the  Poisson  distribution  of  backorders 
at  each  base  for  purchase  alternative  N. 

BPROB 

In  BNEXT,  this  variable  denotes  the  Poisson  probability  that 
exactly  K  spares  will  be  backordered. 

BRATE 

The  expected  number  of  component  removals  per  month  at  each 
base.  Used  only  in  PROLOG. 

A-l 


BRCT 

Base  repair-cycle  time.  The  time  that  elapses  from  the  re¬ 
moval  of  a  base-repairable  unit  to  its  return  to  serviceable 
condition.  Used  only  in  PROLOG. 

BREBO 

In  BNEXT,  the  reduction  in  backorders  that  will  be  experi¬ 
enced  if  an  additional  spare  is  added  to  each  base. 

BRPIPE 

The  base-repair  pipeline.  Used  only  in  PROLOG. 

Cl 

A  constant  used  in  the  subroutine  LIB  in  the  computation  of 
either  depot  or  base  EBO.  It  equals  two-thirds. 

C2 

Another  constant  used  in  the  subroutine  LIB  in  the  compu¬ 
tation  of  depot  or  base  EBO.  It  equals  23/270. 

C3 

Another  constant  used  in  subroutine  LIB  in  the  computation  of 
depot  or  base  EBO.  It  equals  one-sixth. 

CPIPE 

The  condemnation  pipeline.  Used  only  in  PROLOG. 

CRANK(I) 

The  benefit-to-cost  ratio,  or  ranking  value,  of  a  component. 

DBOLIB(N) 

In  MARGIN,  the  value  that  DEBO  will  assume  if  purchase  alter¬ 
native  N  is  best. 

DBOTOT 

DEBO  summed  across  all  components  at  a  particular  avail¬ 
ability  level.  Used  only  in  EPILOG. 

DCOST 

The  total  cost  of  depot  spares  at  a  particular  availability 
level.  Used  only  in  EPILOG. 

DEBO(I) 

The  expected  number  of  backorders  of  a  component  at  the 
depot. 

DELTA 

In  AVAIL,  DELTA  equals  the  difference  between  the  avail¬ 
ability  logarithm  of  the  purchased  component  after  its  pur- 
chase  and  the  availability  logarithm  of  the  component  before 
its  purchase. 

DMEAN(I) 

The  expected  number  of  units  of  a  component  in  shipment  to 
the  depot,  depot  repair,  or  procurement.  This  variable  is 
the  mean  of  the  Poisson  distribution  of  backorders  at  the 
depot. 

DMLOG 

The  log  of  the  depot  pipeline,  DMEAN.  It  is  used  in  the 
computation  of  DEBO. 

DOPIPE 

The  depot  order-and-shipping  pipeline  for  each  base.  Used 
only  in  PROLOG. 

DOST 

The  depot  order-and-shipping  time.  Used  only  in  PROLOG. 

DPIPE 

The  depot  repair  pipeline.  Used  only  in  PROLOG. 
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DPQ 

In  FACTLN,  the  double-precision  equivalent  of  the  variable 
ISIG. 

DPROB(I) 

The  logarithm  of  a  Poisson  probability.  Specifically,  the 
log  of  the  probability  that  the  population  of  the  depot 
pipeline  is  equal  to  JSPARE(I). 

DRBLIB(N) 

In  MARGIN,  the  reduction  in  backorders  that  will  be  provided 
by  an  additional  depot  spare  if  purchase  alternative  N  is 
best . 

DRCT 

Depot  repair-cycle  time.  The  time  that  elapses  from  the 
removal  of  a  depot-repairable  unit  to  its  return  to  service¬ 
able  condition.  Used  only  in  PROLOG. 

DREBO(I) 

The  reduction  in  expected  backorders  that  will  accompany  the 
addition  of  a  depot  spare.  This  variable  strictly  decreases 
as  the  depot  inventory  increases. 

DSPIPE 

The  direct-support  repair  pipeline.  Used  only  in  PROLOG. 

DSRCT 

The  direct-support  repair-cycle  time.  The  time  that  elapses 
from  the  removal  of  a  Direct-Support  repairable  unit  to  its 
return  to  serviceable  condition. 

E 

In  LIB,  CON,  POISSY  or  PRELIM,  the  EBO  that  will  persist  at  a 
site  after  it  has  obtained  a  pipeline's  worth  of  spares. 

GOPIPE 

The  general-support  order-and-shipping  pipeline.  Used  only 
in  PROLOG. 

GSOST 

The  general-support  order-and-shipping  time.  Used  only  in 
PROLOG. 

GSPIPE 

The  general-support  pipeline.  Used  only  in  PROLOG. 

GSRCT 

General-support  repair-cycle  time.  The  time  that  elapses 
from  the  removal  of  a  general-support  repairable  unit  to  its 
return  to  serviceable  condition.  Used  only  in  PROLOG. 

I 

Indicates  a  particular  component.  It  takes  on  values  between 
(and  including)  component  and  NCOMP. 

ICND 

The  percentage  of  parts  condemned.  Used  only  in  PROLOG. 

ICOST(I) 

The  cost  of  a  spare  unit  of  a  component. 

IDONE 

Indicates  whether  all  components  have  been  ranked  in  COMPOP. 
It  is  initially  equal  to  zero;  it  is  set  equal  to  one  once 
the  ranking  process  is  completed. 

IF  IN 

In  POISSY,  the  quantity  of  base  spares  accounted  for  when  the 
iterative  computation  of  EBO  ends. 

IMEAN 

In  LIB  or  CON,  the  truncated  value  of  AMEAN. 
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I  OUT 

When  it  equals  zero,  DRAMA  provides  only  standard  output. 
When  it  equals  one,  the  base  spares  file  is  output.  When  it 
equals  two,  the  depot  spares  file  is  output  as  well.  When  it 
equals  three,  depot  and  base  EBO  files  will  be  output  too. 

IPOS (NCOMP) 

Indicates  which  component  is  in  ranking  position  one,  ranking 
position  two,  etc. 

IPRE 

Dictates  whether  to  call  subroutine  PRELIM.  When  it  equals 
zero,  PRELIM  is  bypassed;  when  it  equals  one,  PRELIM  buys 
conservatively.  When  it  equals  two,  PRELIM  buys  liberally. 

IPS(I) 

The  number  of  parts  installed  per  end  item  or  system. 

IRTS1 

The  percentage  of  parts  repaired  at  the  base  echelon.  Used 
only  in  PROLOG. 

IRTS2 

The  percentage  of  parts  repaired  at  the  direct-support 
echelon.  Used  only  in  PROLOG. 

IRTS3 

The  percentage  of  parts  repaired  at  the  general-support 
echelon.  Used  only  in  PROLOG. 

IRTS4 

The  percentage  of  parts  repaired  at  the  depot  echelon.  Used 
only  in  PROLOG. 

ISIG 

The  truncation  of  the  following  quantity:  The  Base  or  Depot 
pipeline  minus  six  standard  deviations. 

ISP 

In  CON,  the  quantity  of  spares  at  each  base;  a  surrogate  for 
KSP . 

I  STAR 

In  POISSY,  the  quantity  of  base  spares  accounted  for  when  the 
iterative  computation  of  BEBO  begins. 

ISTRT 

In  CON,  the  quantity  of  base  (or  depot)  spares  for  which  an 
EBO  is  computed.  It  equals  either  the  total  inventory  of 
base  (or  depot)  spares  or  the  base  (or  depot)  pipeline  minus 
six  standard  deviations,  whichever  is  less. 

IT 

The  line  in  the  DRAMA  availability-cost  table  that  is  being 
processed  in  EPILOG. 

ITERS 

Counts  the  number  of  marginal  analysis  iterations  in  a  run. 

JJ 

In  DNEXT,  it  indicates  how  many  spares  will  be  redistributed 
in  a  given  purchase  alternative. 

JK 

In  DNEXT,  it  is  the  indicator  of  the  purchase  alternative. 

JSPARE(I) 

The  number  of  spares  of  a  component  at  the  depot. 

K1 

A  constant  used  in  subroutine  FACTLN.  It  equals  one-half  of 
the  square  root  of  TWOPI . 
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K2 

A  constant  used  in  FACTLN.  It  equals  one- twelfth. 

KSP 

In  BOUT,  the  number  of  spares  that  each  base  will  have  after 
a  given  purchase  alternative  is  executed. 

KSPARE(I) 

The  number  of  spares  of  a  component  at  each  base. 

KSTRT 

The  number  of  base  spares  that  BNEXT  has  accounted  for  in  the 
process  of  computing  BBOLIB. 

K(NSUBS) 

The  subsystem  "K-f actor."  The  ratio  of  subsystem  removals  in 
the  DRAMA  data  base  to  the  subsystem  removals  in  the  raw  data 
base.  Used  only  in  PROLOG. 

LCOST 

The  cost  of  a  particular  component's  spares  in  EPILOG. 

LLL 

In  PRELIM,  the  truncated  value  of  PM. 

LNFACT (30) 

In  FACTLN,  an  array  of  the  logarithms  of  the  factorials  of  1 
through  30. 

M 

In  BOUT,  the  indicator  of  the  purchase  alternative  being 
considered. 

MAXRED 

The  maximum  number  of  parts  that  may  be  redistributed  from 
the  depot  to  the  bases  in  a  given  purchase. 

MODE 

When  it  equals  one,  MARGIN  computes  a  backorder-reduction  to 
cost-ratio,  instead  of  a  marginal  availability-to-cost  ratio, 
for  each  purchase  alternative.  It  is  set  equal  to  one  only 
when  a  component's  marginal  availability  is  undefined,  e.g., 
when  it  has  more  backorders  at  each  base  than  applications. 

N 

In  MARGIN,  it  indicates  the  purchase  alternative.  It  takes 
on  values  from  1  to  NOB+1. 

NAVAIL 

The  line  in  the  DRAMA  availability-cost  table  that  is  to  be 
output  by  EPILOG.  When  it  is  equal  to  zero,  all  lines  are  to 
be  output. 

NBEST(I) 

The  purchase  alternative  which  will  yield  the  most  avail¬ 
ability  for  its  cost. 

NCOMP 

The  number  of  different  types  of  line-replaceable  units 
(LRUs)  in  an  end-item. 

NEWPOS 

In  COMPOP  the  new  ranking  position  of  the  component  that  was 
just  purchased. 

NOB 

The  number  of  bases. 

NT ABLE 

The  number  of  cost-availability  pairs  that  have  to  be  input 
to  EPILOG. 
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p 

In  LIB,  the  Poisson  probability  that  the  quantity  of  parts  in 
the  depot  (or  base)  pipeline  is  exactly  equal  to  the  mean  of 
the  pipeline. 

PLT 

Procurement  lead  time.  The  time  that  elapses  from  the  con¬ 
demnation  of  a  component  to  the  receipt  of  a  replacement  at 
the  depot. 

PO 

In  PRELIM,  LIB,  or  CON,  the  temporary  value  of  BMLOG  or 
DMLOG . 

Q 

In  LIB,  an  intermediate  value  in  the  computation  of  depot  or 
base  EBO. 

QRANK(N) 

In  SIZOP,  the  marginal  availability  (or  marginal  EBO  reduc- 
tion)-to-cost  ratio  for  purchase  alternative  N. 

R 

In  LIB,  or  CON,  the  reduction  in  BEBO  that  will  accompany  the 
addition  of  one  spare  to  each  base. 

RATE 

The  frequency  of  component  removals.  The  number  of  removals 
per  component  application  per  hour  of  use.  Used  only  in 
PROLOG. 

SPB 

The  number  of  systems  (or  end  items)  per  base. 

STIRL 

In  FACTLN,  LIB,  and  CON,  the  logarithm  of  the  factorial  of 
KSP ,  where  KSP  is  the  number  of  spares  accounted  for  at  each 
base . 

SYSLOG 

The  logarithm  of  the  system's  availability. 

SYSRAT 

The  system  utilization  rate.  Expressed  in  usage  hours  per 
month.  Used  only  in  PROLOG. 

SYSTAV 

The  system  availability. 

TAIL 

In  LIB,  the  area  under  the  Poisson  probability  density  func¬ 
tion  between  IMEAN  and  infinity. 

TCOST 

The  total  cost  of  spares  at  a  particular  availability  level. 
Used  only  in  EPILOG. 

TEMP(NCOMP) 

In  COMPOP,  the  component  that  is  temporarily  ranked  number 
one,  number  two,  etc. 

TOTCST 

The  total  cost  of  spares  of  all  component  types. 

TWOPI 

Another  constant  in  the  subroutine  LIB.  It  equals  two  times 
Pi. 

Y 

In  LIB,  an  intermediate  value  in  the  computation  of  base  or 
depot  EBO. 
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APPENDIX  B 
PROGRAM  LISTINGS 


PROLOG 
1  ist 

1000  REAL  K ( 30 ) 

1010  CALL  FMEDIA<9,0> 

.1 020  CALL  FMED I A  <  1 0 ,  0 ) 

1 030 C 

1  040C*****SCENAR  1 0  DAT  A***** 

1 050 C 

1 060  READ (9,1 09 )  NOB , SPB, NSUBS , NCOMP , SYSRAT 

10700 

1 030C*****SUBSY3TEM  DAT A***** 

1 0900 

1100  DO  1000  ISS  =  1, NSUBS 

1110  1 000  READ ( 9 ,  1 09 )  K  < I S3 ) 

1120  109  FORMAT (V ) 

1  1300 

1 140C*****ASSIGN  VALUES  TO  THE  REPAIR  AND  SHIPPING  TIMES.***** 


1  1 50  C 

1 1 60C 

***NORMALLY ,  THESE  TIMES  ARE  EXPRESSED  IN  MONTHS*** 

1  1 70  C 
1  130C 

**NOTE: TO 

MODEL 

A  PERFECT  TWO-ECHELON  SYSTEM,  SET 

1 1 90  C 

OSRCT,  G: 

3RCT , DS 

OST,AND  GSOST  TO  ZERO. 

1 200C 
1210 

BRCT  =  . 1 

1  ■  1  t  ) 

OSRCT  =  .5 

1 230 

GSRCT  =  1.0 

1240 

DRCT  =  5.0 

1230 

PLT  =  18.0 

1260 

GSOST  =  . 5 

1 270 

DOST  =  2. 0 

1 23  OC 
1 2900 

****«WR I TE  SCENAR I 0 

INPUT' 

TO  DRAMA  INTO  FILE  CODE  10.***** 

1 3000 

1 3 1 0  WR I TE ( 1 0 , 1 09 )  NOB , SPB , NCOMP 

1 3200: 

1 330C*****N0W  INPUT  THE  COMPONENT  DATA.  PROCESS  IT.  THEN  OUTPUT 
13400  I  F*S ,  I  COST ,  BMEANO ,  AND  DMEAN  TO  FILE  CODE  10. 

1 35  OC 

1 360  DO  2000  1=1, NCOMP 

1 370  READ ( 9 ,  1 09 )  I COST , I PS ,  I RTS 1 ,  I RTS2 ,  I RTSS , I RTS4 -  I CND , I  S3 , RATE 

1 38  OC 

1 3'~'0C  *-****COMPUTE  BRATE  FOR  THIS  COMPONENT .  ***** 

1400  BRATE  =  K(ISS)  *  RATE  *  SYSRAT  *  FLOAT (IPS)  *  SPB 

1 4 1  OC 

1  4200  * * * * * C 0 M P U T E  DMEAN  FOR  THIS  COMPONENT.  ***** 

1430  DP  I PE  =  FLOAT (NOB)  *  FLOAT  < I RTS4/ 1 00 )  *  BRATE  *  DRCT 

1440  CP  I PE  =  FLOAT (NOB)  *  FLOAT ( I CND/ 1 00 )  *  BRATE  *  PLT 

1 430  GSPIPE  =  FLOAT ( NOB )  *  FLOAT ( I RTSS/ 100)  *  BRATE  *  GSRCT 

1460  DMEAN  =  DP  I PE  +  CP  I PE  +  GSPIPE 


1  4  7  OC 

1  4  SOC: 

* *** *COMPUTE  BMEANO  FOR  THI 

]Z; 

C  0  M  P  0 N  ENT. *  * * * * 

1  490 

BRP I PE  =  FLOAT ( I RTS  1 / 1 00 ) 

BRATE  *  BRCT 

1 300 

DSP  I PE  =  FLOAT ( I RT 32/ 1 00 ) 

<- 

BRATE  *  OSRCT 

1  3  l  0 

GOP  I PE  =  FLO A T ( I RTS 3 / 1 00 ) 

-X- 

BRATE  *  GSOST 

B-l 


1520  OOP I PE  =  FLOAT <  < IRTS4+ICND) / 100  )  *  BRATE  *  DOST 

1530  BMEANO  =  DMEAN/ FLOAT  (  NOB  )  +  BRPIPE  +  DSP  I  PE  +  GOP  I  PE  +• 

1540  &  DOF' I  PE 

1 5500 

1560C**hhh»0UTPUT  COMPONENT  DATA,  *■*■**■* 

1 570  2000  WR I TE  < 1 0 , 1 09 )  I PS , I COST , BMEANO , DMEAN 
1580  STOP 

1590  END 
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DRAMA 


LIST 


i  000 
1 0 1 0 
1 020C 
1 030 
1  Q&Q 
1  050 
1  060 
1070 

1  oso 

1 090 
1  1  00 
1110 
1  120 
1  1 30 
1 1 40 
1  1 50 
1  1 60 
1  1 70 
1  1 80C 
1 1 90  C 

1 20  OC 

1 2 1  OC 
1  220 
1 230 C 
1 240C 
1250t 
1 2  60  ( 

.1 2  TOC 
1 2S0  C 
1 29  OC 
1 300i 

1 31  OC 

1 32  OC 
1 330 
1340 
1  35  OC 
1 36  OC 
1 3  TOC 
1  380C 
1 3  9  Of 


COMMON/ AVA/ AVLOG C 70  > , CRANK ( 70 ) ,  IPS < 70 ) 

Ci  MM  ON  /  BASE  /  BEBO 

COMMON/EBON/BEBON ( 70 ) , DEBON <  70 ) 
OMMON/CONST/Cl * C2, 03, TWOPI 
OMMON /DEPOT /DEE'CK  DMLOO,  DPROB,  DREBO 
0MMGN/KQNST/K1 * K2 

OMMQN/NEXT/BBOLIB, OBOLIB, PROLIB? DRBLIB 
OMMON /MEAN/BMEANO ( 70 ) , DMEAN ( 70 ) 

OMMON /PAR AMS/ I DONE  ? MAX RED , MODE, NCOMP? NOB, SPB 
OMMON /PR  I CE/ 1 COST ( 70 ) , NBEST C 70 ) 

PARE C  70 ) , K SPARE C  70 > 

LOGS  SYSTA'vS  TOTCST 


OMMON/ SPARE /JSf 
IS  M  M  0  M  /  3  Y  S  T  /  S  Y  S  l 
OMMON /OPTS/ I PRE 
REAL  BEBO (70) 

REAL  DEBO  ( 70 )  ,  DMLOG  (  70 )  ,  DPROB  (  70  > 

DOUBLE  PR EC I S I ON  C 1 , 02 , 03 , TWOP I 
DOUBLE  PR EC I S I ON  K 1 , K2 

DOUBLE  PR  EC  1 3 1  ON  BBOL  IB(  13)  ,  DBOLI B  (  1 3  )  ,  PR  01. 1 E  (  1 3  )  ,  DRBL I B  (  1 3  ) 


*****THIS  SECTION  OF  DRAMA  ATTACHES  THE  NECESSARY  DATA  FILES. 


DREBO ( 70 ) 


THE  FIRST 


.E  ATTACHED  IS  THE  INPUT  FILE. 


'HE  LAST  FILES 


T ACHED  ARE  OPTIONAL  OUTPUT  FILES. 

ALL  ATTACH  <  9 ,  **  0S29  /  N232D  /  LEMON  ?  ”  ,  3 ,  0 ,  I  ST  AT ,  ) 

ALL  FMEDIA< 10,0) 

ALL  FMED I A ( 1 1 »  0  > 

ALL  FMED I A < 12,0) 

ALL  FMED I  ft < 13, 0) 


***THI3  SECTION  OF  DRAMA  ASSIGNS  INITIAL  VALUES  TO  ALL  KEY 
KEY  VARIABLES.  USERS  MUST  ASSIGN  VALUES  HERE  IF  THEY 
WISH  TO  EXECUTE  OFT  IONS 

*****USER  DEFINED  CONTROL  VARIABLES***** 

I  PRE  =  0 
I OUT  =  0 

*****NGTE:  IF  THE  USER  SETS  I OUT  TO  3  OR  4  HE  MUST  DO  TWO 

THINGS  IN  ORDER  TO  OUTPUT  EXACT  VALUES  OF  BEBO  AND  DEBO. 
ONE ,  HE  MUST  DEFINE  "BEBON"  AND  "DEBON"  IN  THE  COMMON 
SECTIONS  OF  DRAMA  (SEE  LINE  1020),  MARGIN,  AND  UPDATE. 
TWO,  HE  MUST  RECOMPUTE  BEBO  IN  SUBROUTINE  UPDATE  INSTEAD 
JBROUTINE  MARGIN.  SEE  COMMENTS  BELOW. 

OQQ 


1 400 C 

OF 

1410 

AVLIM 

1420 

COLIM  : 

1  OC: 

1  44.00 

«•####  I N 

1450 

ITERS 

1  460 

I  DONE 

1470 

AVTEST 

1  4300 

1  4900 

-K-  -&• Y 

1 500 

READ ( 9 

1 5 1  OC 

C 
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.1  52  OC 
1  530 

1540  900 
1350  109 
1  56  OC 
1  57  OC 
1  58  OC 
1 59  OC 
1  600 
1  6 1 0 
1  620 
1 630 
1  640 
1 650 
1 660C 


<- COMPONENT-  LEVEL  VAR  I  ABLES***** _ 

DO  900  1=1, NCOMP 

READ ( 9 , 1 09 )  I PS ( J ) , 1 COST ( I ) , BMEANO ( I ) , DMEAN ( I ) 

FORMAT ( V ) 

***NOTE:  IN  THIS  FORMAT,  ONE  SPACE  SEPARATES  EACH  DATA 
ELEMENT  FROM  ITS  SUCCEEDING  ELEMENT. 


** ***SUBROUT  I NE  CONSTANTS***** 
f:  i  =  2  -  0  Ei  0  /  3 . 0  D  0 
C  2  =  2  3 ..  0  D 0  /  2 7  0 . 0  D  0 
Co  =  . 25D0  *  Cl 

TWOP 1=2- ODO  *  3 .141 592 6 5358 97 32D0 
K1  =  .500  *  DLOG(TWGPI) 

K  2  =  .  5D0  ->  C3 


1670C**#**WRITE  OUTPUT  HEADER 
1630  WRITE (6, 996 ) 

1 690  996  FORMAT  (  ***  '*  ,  '  A V A I  LAB  I L I  TV  "V  1  6  X  ,  *'  COST  "*  ) 

1 700  C 
1 7 1  OC 


1  72  OC 
1 730C 
1 740C 
1  750C 
1 7 60 C 
1  770 
17S0C 
1 790C 
1 300C 
1 8 1  OC 
1  320C 
1 330 C 
1 340 

1 35  OC 

1 36  OC 


*****THI S  SECTION  OF  DRAMA  INITIALIZES  THE  VARIABLES  INVOLVED  IN 
THE  COMPUTATION  OF  PESO  AND  DEED.  THE  DO  LOOP  BELOW  ALSO 
INITIALIZES  COMPONENT  AVAILABILITIES  AND  COMPUTES  THE  INI¬ 
TIAL  COMPONENT  RANKING  VALUES. 


DO  1000  I  =  1, NCOMP 

***NOTE :  THE  VARIABLE  "MODE"  TELLS  THE  SUBROUTINE  MARGIN 
WHETHER  A  COMPONENT  HAS  A  POSITIVE  AVAILABILITY.  IF  fi 
COMPONENT  HAS  A  NON-POSITIVE  AVAILABILITY,  THEN 
MODE  STAYS  EQUAL  TO  2,  AND  AN  AUTOMATIC  PRELIMINARY 
PURCHASE  WILL  BE  MADE  IN  ORDER  TO  MAKE  AVAILABILITY- 
POSITIVE. 

MODE  =  1 


***** INITIALIZE  SPARES  I N VENTOR I ES***** 


1  370 

JSPARE(I)  =  0 

1 330 

KSPARE(I)  «  0 

1 39  OC 

1 900 C 

***** INITIALIZE  THE  VARIABLES  INVOL1 

VED  IN  THE 

EBO  COMF'UT A— 

1 9 1  OC 

TIONS 

1 920 

BEBO(I)  =  BMEANO ( I ) 

1 930 

DEBO(I)  =  DMEAN ( I ) 

1  940 

DMLOG < I )  =  DLOG <  DMEAN ( I ) ) 

1 950 

DF'ROE: (  I )  =  tDMEAN(I) 

1  960 

DREBO(I)  =  l.ODO 

1  970C 

***NOTE :  THE  INDEX  OF  A  DO-LOOP 

CANNOT  BE 

PASSED  TO 

1 930C 

A  SUBROUTINE.  THEREFORE,  IMAX 

MUST  BE  U: 

SED  AS  A  SUR- 

1 990C 

GATE  FOR  " I " . 

2000 

IMAX  =  I 

20100 
20200 
2030 
20400 

2050C*****  INITIALIZE  TH I  S  COMPONENT S  AVA I  LAB  I L I  TY***«* 


***MAKE  A  PRELIMINARY  BUY  OF  THIS  COMPONENT  IF  ONE  IS  CALLED  FC 
IF  <  IPRE  . GT.  0  .AND.  DMEAN ( I )  . GE.  10.0  )  CALL  PRELIM ( IMA 
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CALL  AVAIL (I MAX) 


>060 
1:070c 
I'OSOC 
•090i 
>  1 00 
• 1  1 OC 
1 1 20C 


*****COMPUTE  THIS  COMPONENT- S  RANKING  VALUE***** 
1 000  CALL  MARGIN ( I MAX ) 


2 1  30  C 

*****GOMPLETE  THE  FIRS 

T  ITERATION  i 

2 1  40C 

THE  COMPONENTS  AND  IDENTIFY  T! 

2 1  50C 

" I DONE" . 

2 1 60 

CALL  CO MR OP ( I MAX 

) 

2170 

I DONE  =  1 

2 1 SOC 

2 1 90  C 

*****PURCHASE  THE  BEST 

COMPONENT  **■ 

2200 

2000  CALL  UPDATE (IMAX) 

2210 

ITERS  =  ITERS  + 

1 

2220C 

*****TEST ; SHOULD 

AVAILABILITY 

223 OC 

THIS  POINT? 

2240 

IF  (  SYSTAV  .LT. 

AVTEST  )  00 

2250 

IF  (  AVTEST  . OE. 

. 989  )  AVTE 

2260 

IF  <  AVTEST  .OE 

. 89  . AND .  AV 

— .L  /  O 

Zt  AVTEST 

+  ,  0 1 

2230 

IF  (  AVTEST  .OE. 

•  V •_*  »  i  ri 

2290 

?■:  AVTEST  *  .0 

cr 

2300C 

*****OUTPUT  DATA 

IF  THE  TEST 

2  6 1 0 

WRITE! 6, 1006)  SY 

STAV , TOTCST 

2320 

1 006  FORMAT  (  ,  F 1 2 .  S 

, 1 0  X , F10.0) 

2330C 

*****OPT I CNAL  OUT 

PUTS***** 

2340 C 

IF  (  I out  :gt. 

( j  )  HR I TE ( 1 0 

23  SOC 

IF  (  TOUT  .07. 

1  !  WR I TE  < 1 1 

2360C 

IF  (  I OUT  . OT. 

2  )  WRITE <12 

2370C 

IF  (  I OUT  .OT. 

2  i  WRITE <13 

2330C 

2390C 

2400C 

***** TEST :  HAS  A  CONST 

■RAINT  BEEN  R 

2410 

3000  CONTINUE 

2420 

IF  <  SYSTAV  . OE. 

AVLIM  .OR. 

.  LT. 


)  A V TEST  = 


2430 

,1  .•'!  r 

_ -r  ~T  V  «. 

L  4  5  O 

246*0 

2:470 

2480C 

2490 

2500 

2510 

252 OC 

2530 

2540 

2550 

2560C 

2570 

2  5  SO 

2590 


O  CO  70  -6000 

:**«**  IF  NEITHER  CONSTRAINT  HAS  BEEN  REACHED  REPEAT  MARC 
CALL  MARGIN (I MAX) 

CALL  CO MR OF ( IMAX ) 

GO  T 0  2000 

:***** IF  EITHER  CONSTRAINT  -HAS  BEEN  REACHED,  STOP.***** 
6000  CONTINUE 
STOP 
END 

SUBROUT  I NE  MARC I N ( I ) 


COMMON/ 


AVA/  AVLril';  r  7 


00 ( 70 ) , CRANK ( 70 ) , I PS ( 70 ) 


COMMON /BASE/ BEBO 
COMMON / EBON / BEBON ( 70 ) , DEBON ( 70 ) 

COMMON /DEPOT  /DEBO ,  DM  LOG’ ,  DPROB ,  DREBO 
COMMON/NEXT/BBOL  I  B ,  DBOL  I B -  F'ROL  I B ,  DRBLI  B 
COMMON/PARAMS/ I  DONE , MAX RED , MODE - NCOMP ,  NOB ,  S 
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ALL 

i^ET 


1*600 

1610 

1620 


2640< 


>C 


26601 
2670 C 
26  SO 
26 900 
2700C 
27 1 00 
2720  2010 
2730 
27400 
2 7500 


i 


27700 
27S00 
2790 
280C  “ 
28 1  < 
28201 
2830C 
2840C 

23 6< 
2870 
28801 
2890 
2900 
29 1 OC 
2920 C 
2930C 
2940C 

295< 
2960 
2970 
2980 
2990 
3000 
30 1 0 
3020 


COMMON/PR I CE/ I COST ( 70 ) , NBEST  <70 ) 

REAL  BE BO  <  70 ) 

REAL  DEBQ ( 70 ) , DMLOG ( 70 ) , DPROB < 70 ) , BREED <  70 ) 

DOUBLE  PREC I S I ON  BBOL I B  < 1 8 ) , DBOL I B ( 1 8 ) , PROL I B (1 8 ) , DRBL I B  < 1 8 ) 


;*****TEST:  CAN  BEBO  BE  FURTHER  REDUCED?  .  IF  IT  CANNOT, 

SKIP  BEBO  CONFUTATIONS  AND  SET  THE  RANKING  VALUE  TO  ZERO. 


IF  (  BEBO < I )  .  LE.  0.0  )  GO  TO  1510 
##**-*COliPUTE  BBOL  I B^*-*** 


CALL  DNEXT < I ) 
CALL  BNEXT(I) 


;*****BEBG  ■  OUTCOMES  ARE  NOW  STORED  IN  THE  ARRAY  "BBOL IB”.  USE 
THEN  TO  CONFUTE  THE  BENEFIT-TO-COST  RATIOS  “OF  THE  FEAS- 
::  I BLE  PURCHASES .  THEN  FIND  THE  MOST  COST-EFFECTIVE  ONE. 


CALL  SIZOP(I) 


;*****THE  MOST  COST  EFFECTIVE  PURCHASE  HAS  BEEN  DESIGNATED 
" NBEST" .  STORE  THE  DEBQ  AND  BEBO  VALUES  ASSOCIATED 
::  WITH  '=  :S  PURCHASE. 


BEBO(I)  =  BBOLIB<  NBEST < I )  ) 
BEBON < I )  =  BBOL I B ( NBEST  < I ) ) 
DEBO ( I )  =  DBOL I B (  NBEST  < I )  ) 

DEBON < I )  =  DBOL I B ( NBEST ( I ) ) 
DPROB ( I  )  =  PROL I B  <  NBEST ( I  )  ) 
DREBO(I)  =  DRBL IB <  NBEST ( I  )  ) 


::**•■*-** IF  THIS  COMPONENT  ALREADY  HAS  A  DEFINED  AVAILABILITY,  RETURN. 
IF  IT  DOES  NOT,  BUY  SPARES  UNTIL  ITS  AVAILABILITY  CAN  BE 
COMPUTED. 


IF  <  MODE  .  EGi.  2  )  GO  TO  1610 
CALL  UPDATE  < I ) 

GO  TO  2010 

1510  CRANK (I)  =  0.0 
1610  CONTINUE 
RETURN 
END 


3040  SUBROUTINE  DNEXT ( I ) 

3050  COMMON/NEXT/BBOLIB, DBOLI B , PROLIB, DRBL IB 

3060  COMMON / DEPOT / DEBO  > DMLOG, DPROB, DREBO 

3070  COMMON / SF ARE/USF ARE ( 70 ) , KSPARE <70 ) 

3080  COMMON/ P ARAMS / I DONE , MAXRED, MODE, NCOMP, NOB, SPB 

3090  REAL  DEBO < 70 ) , DMLOG  < 70 ) , DREBO <  70 ) , DPROB <  70 ) 

3 1 00  DOUBLE  PREC I S I ON  BBOL I B  < 1 8 ) , DBOL I B  < 1 8 ) , DRBL I B  < 1 8 ) , FROL I B  < 1 8 ) 

31  IOC 

3120C*****C0MPUTE  THE  DEBO  THAT  WILL  RESULT  IF  THE  DEPOT  GETS  A 
3130C  SPARE.  THIS  IS  PURCHASE  ALTERNATIVE  NOB+1. 
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cn  ui  '  n  cn  i 


3 1  40  C 
3 1 30 
3 1 60 
3170 
3 1  30 
3 1 90 
32000 

3210C*****ST0RE  THE  CURRENT  VALUES-  OF  DE30,  DPROB,  AND  DREBO 


N  =  NOB  +  1 

DRBL I B  (  N )  =  DREBO  (  I  )  -  DEXP  (  DF'ROB  (  I  )  ) 

DBOLIB ( N )  =  DEEO(I)  -  DRBLIB (N) 

FROL I B ( N )  =  DPROB ( I  )  -  DLOO (  1 . ODO  +  DBLE ( JSP ARE < I )  )  )  + 

DMLOG < I ) 


THIS 


322 OC 

PURCHASE  ALTERNATIVE  NOB. 

3230  C 

3240 

DBOL I B ( NOB ) 

—  DEBO ( I ) 

3250 

DRBL IB (NOB) 

=  DREBO < I ) 

3260 

PROL I B ( NOB ) 

=  DPROB ( I ) 

327 OC 

3280 C 

***H**CONPUTC  THE  VALUES  OF  DEBCu  DPROB 

329 OC 

UTHER  PURCHASE 

ALTERNATIVES. 

LL  OF  THE 


DO 


***** DETERM  I NE  THE  MAXIMUM  NUMBER  OF  SPARES  THAT  CAN  BE  RE¬ 
DISTRIBUTED.  IF  NONE  CAN  BE  REDISTRIBUTED,  RETURN. 


MAXRED  =  M I NO (  NOB  -  1,  JSPARE(I)  ) 

IF  (  MAXRED  .  EQ.  0  )  00  TO  5015 

i-****PERFORM  COMPUTATIONS  FOR  EACH  PURCHASE  ALTERNATIVE- 
ALTERNATIVE  NOB-1  FIRST,  THEN  DO  NOB-2,  NOB-3,  ETC. 

DO  4415  JJ  3  1,  MAXRED 
J  =  NOB  -  JJ 
JK  =  J  +  1 

DBOL I B  <  J )  =  DBOLIB ( JK )  +  DRBL I B  C JK ) 

FROL I B ( J )  =  PROL I B ( JK )  +  DLOO ( DELE ( JSP ARE ( I ) )  - 
DBLE ( JJ )  +  1 . ODO  )  -  DMLOG ( I ) 

DRBL  IB  ( J  )  =  DRBL  I B  (  JK )  4-  DEXP  ( PROL  I B  ( J  )  ) 

I F  (  DRBL I B ( J )  . OE .  1 . ODO  )  DRBL I B ( J )  =  1 . ODO 


33 00 C 
33 1 OC 
3320C 
3330 C 
3340 
3350 
3360 C 
3370C 
3  3S0C 
3 39 OC 
3400 
3410 
3420 
3430 
3440 
3450 
34  ao 

3470  4415 
3430 
3490 
3500 
35 1  OC 
3520 
352:0 
3 5 40 
:":F50 
60 
70 
SO 
90C 
3600C 
36 1  OC 
3620 
3630C 
36  4  OC 

3650C*****CGMPUTE  THE  BASE  EXPECTED  BACKORDER  (BBOLIB)  VALUE  FOR 
3660C  EACH  PURCHASE  ALTERNATIVE. 

■3670 C 


5015  CONTINUE 
RETURN 
END 


SUBROUTINE  BNEXT(I) 

COMMON/NEXT/BBOLIB, DBOLIB, PROL IB, DRBLIB 
COMMON/MEAN/ BMEANO  < 70 ) , DMEAN < 70 ) 

COMMON /SPARE / JSP ARE ( 70  > , KSPARE ( 70 ) 

COMMON /PAR  AMS/  I  DONE-,  MAXRED,  MODE,  NCOMFs  NOB,  SPB 
DOUBLE  PREC I S I ON  BEOL I B  < 1 S ) , DBOL I B (1 8 ) , PROL I B ( 1 3 )  ,  DRBL I B ( 1 S ) 
DOUBLE  PREC I S I ON  BMEANN , BPROB , BREED , BMLOG 

*****COMPUTE  THE  BASE  SPARES  INVENTORY  FOR  THE  NEXT  PURCHASE 

KSP  =  KSPARE  (I)  4-  1 
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i 


36S0 

3690 

3700C 

371  OC 

372 OC 

3730 

3740 


3800 


33 1  OC 
3 32 OC 
333 OC 
3  340  C 
3S50C 
33-’tOC 
337 OC 
383 OC 
3590  C 
3900 
39 1  0 
3920 
3930 
3940C 
3950C 

3  9  tOC 
3 97  OC 
393  OC 
399 OC 
400 OC 
40 1  0 
4020 
403 OC 
4040C 
4050C 
4060 C 
4 0 7 0 C 
40S0C 
4090 
4 1 00 
411  OC 

4  1 20C 
4  1 30C 
4140C 
4  1 50C 
4  1  600 
4  1  7  OC 
4 1  SO 
4 1 90 
4200  C 


DC  2020  N  =  NOB  -  MAXRED*  NOB  -  1 

IF  (  N  . EQ.  NOB  +  1  )  KSP  =  KSP  - 

*****RE- INITIALIZE  BMEAN  FOR  EACH  PURCHASE  OUTCOME. 

BMEANN  =  DBLE ( BMEANO ( I ) )  -  (  DBLE ( DMEAN ( I ) >  - 
DBOL I B ( N )  ) / DBLE ( NOB ) 

IF  THE  BASE  HAS  NO  SPARES »  AND  IF  IT  IS  ALLOTTED  NO 
SPARES  IN  THIS  OUTCOME,  THEN  BBOLIB  EQUALS  BMEAN. 

BBOLIB <N)  =  BMEANN 

IF(  KSP  .El'!.  0  )  GO  TO  2020 


♦♦♦♦♦IF  THE  BASE  HAS  SPARES,  BBOLIB  WILL  BE  LESS  THAN  BMEAN. 
COMPUTE  BBOLIB  FROM  BMEAN  GIVEN  KSP  BASE  SPARES. 

DO  THIS  IN  THREE  STEPS.  FIRST,  INITIALIZE  KEY  VARIABLES. 
SECOND,  USE  SHORTCUT  CALCULATIONS,  IF  POSSIBLE. 

THIRD,  USE  AN  ITERATIVE  COMPUTATION,  IF  NECESSARY. 

*****INITIALI ZE 

KSTPT  =  0 

BMLOO  =  DLOG( BMEANN) 

BP ROB  =  -BMEANN 
BREBO  =  1 . ODO 

-s-  *•**•«■  SHORT  CUTS***** 


SHORTCUT  1  COMPUTES  THE  VALUE  OF  BBOLIB  FOR  AN  IN¬ 
VENTORY  OF  BMEAN.  THEN  IT  SETS  KSTRT  EQUAL  TO  THIS 
I  INVENTORY  QUAN  T I TY  . 

IF  <  KSP  . GE.  IFIX( BMEANN 5  .AND.  BMEANN  . GE.  10.0) 

CALL  L I B ( BP ROB , BREED , BBOL I B < N ) , BMLOO , BMEANN , KSTF 

♦♦SHORTCUT  2  COMPUTES  THE  VALUE  OF  BBOLIB  FOR  AN 
INVENTORY  OF  KSP,  OR  BMEAN  MINUS  SIX  STANDARD  DE¬ 
VIATIONS,  WHICHEVER  IS  LESS,  THEN  IT  SETS  KSTRT 
EQUAL  TO  THIS  INVENTORY  QUANTITY 

IF  <  KSP  . LT.  I  FIX (BMEANN)  .AND.  BMEANN  . GE.  SS.O) 
CALL  CON ( BP ROB , BBOL I  3 ( N ) , BMLOO , BMEANN , KSTRT ,  KSP ) 


♦♦♦♦♦ ITERATIVE  COMPUTATION  OF  BBOLIB. 


IF  THE  INVENTORY  IS  GREATER  THAN  KSTRT,  ONLY  PART  OF 

THE  INVENTORY  HAS  BEEN  ACCOUNTED  FOR.  THUS 

POISSY  MUST  BE  CALLED  TO  ACCOUNT  FOR  THE  REMAINDER. 


IF 


(  KSP  .GT.  KSTRT  )  CALL  POISSY  (  BF'ROB,  BREBO, 
BBOLIB ( N )  ,  BMLOO , KSTRT , KSP ) 


0  CONTINUE 
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RETURN 

END 


4220 
4230 
4240C 

4230  SUBROUT  I ME  L I B ( P , R , E , FO , AM , I MEAN ) 

4260  COMMON/ CONST /C 1 , C2 , C3 , TWOP I 

4270  DOUBLE  PREC 1 3 1  ON  C 1 , C2 > C3 , TWOP I 

4230  DOUBLE  PREC 1 3 1  ON  P , R , E , PO , AM , Y , Q .TAIL, ST  I RL 

4290C 

4300  I MEAN  =  I FIX (AM) 

43 1 OC 

4320C*****C0MPUTE  THE  LN  OF  1MEAN  FACTORIAL***** 

4330C 

4340  CALL  FACTLN ( 3TIRL , IMEAN ) 

4 3300 

4 3 6 0 C***** COMPUTE  THE  POISSON  DPROBABILITY  THAT  IMEAN  UNITS  ARE  IN  THE 
437 OC  PIPELINE. 

4330C 

4390  P  =  DBLEt IMEAN)  *  PO  -  AM  -  STIRL 

4400C 

44 1 OC*****COMPUTE  7HE  AREA  UNDER  THE  PDF  FROM  IMEAN  TO  INFINITY. 

4420C  USE  THE  ALGORITHM  FROM  KNUTH,  VOLUME  1.  F.  116. 

443 Oi  C 


AM  -  DELE (IMEAN) 

Y  —  Cl  +  (  C2  —  C3  *  Y  /  2.0 DO  —  C3  *  Y  *  Y  *  Y 
DELE ( IMEAN) 

=  Q  /  DSQRT (  DBLE( IMEAN)  *  TWOP I  )  +  . 3D0 
AREA  UNDER  THE 


4440 

\  /  _ 

T  — 

Q  s 

4460 

4470 

TAIL 

448  OC 

4490C 

»****R  IS  THE 

4300 C 

CLUB  I  NO 

45  IOC 

4320 

R  = 

4530C 

45400 

*****COMPUTE 

4530C 

4560 

E  =  Y  * 

4570C 

4330 

RETURN 

4590 

END 

THE  PD  OF  IMEAN. 
TAIL  +  DEXP(P) 
THE  EEC  AT  THIS 
TAIL  +  AM  *  (  R 


PDF  FROM  IMEAN 

SITE  GIVEN  IMEAN 
-  TAIL  ) 


TO  INFINITY  IN¬ 


SPARES  . 


)  / 


4600C - 

46 1  0  SUBROUT I NE  CON  <  P , E , PO , AM ,  I STRT , I SP ) 

4620  DOUBLE  PREC 1 3 1  ON  P , E , PO -  AM , ST I RL 

4630)  C 

4 6 4 0 C  *  *  *  *  *  C  0  M  P  U  T  E  THE  'INTEGER  THAT  IS  SIX  STD.  DEVIATIONS  BELOW  THE  MEAN. 
46500 

4660  1SIG  =  I F I X  <  AM  -  6.0  *  SORT (AM)  ) 

467  OC 

46S0C  *****SELECT  THE  SPARES  QUANTITY  FOR  WHICH  EBO  WILL 

4690C  BE  CALCULATED. 

4700C 

4710  I STRT  =  MI  NO  (  ISIG,  ISP  ) 

4720C 

4730C*****DETERMINE  UNSATISFIED  LI  BAND. 

4740)0 

4730  E  =  AM  -  DELE (I STRT) 
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4760 

4770C 

IF  <  ISP  . LT. 

4730C 

***** IF  THE  SPARES 

4790C 

□PROBABILITY 

4300C 
43 1  DC 

TO  ISIG. 

432 OC 
433 OC 

*****COMPUTE 

4  Q40 
4350C 

CALL  FAC 

436  OC 
4370 C 

****  *COMPUTE 

COMPUTE  THE  POT 


4330 

p  = 

DBLEC IS 

IG)  *  pi 

439 OC 

4900 

302: 0 

CONTINUE 

4910 

RETURN 

4920 

END 

493 OC 

— 

4940 

SUBROUT I 

NE 

FACTLN  <  ST  I  R| 

4950 C 

-X.  -X.  J<.  -X.  -X.  -X 

■*  THIS  S 

IJBROUTI 

NE  GENE. 

496  OC 

-X-  -£•  <- 

*  IF  N  I 

<*• 

LESS 

THAN  31 

497 OC 

STORED 

I 

N  THE 

TABLE  , 

493  OC 

4990 

COMMON /K 

ON 

ST/K1 

K 

5000 

DOUBLE  P 

RE 

Cl  SION  DPQ, S 

50  1 0 

DOUBLE ■ P 

RE 

C I S  1 0 

N  LNF  AC 

5020 

LiATA  LNF 

AC 

7/ 

503 0 

*. 

0 . ODO » 

5040 

.6901471 

so 

55994 

53 1 ODO ^ 

5050 

?/ 

. 1791759 

46 

92230 

5500 D  1 , 

5060 

. 3 1 78053 

03479 

4562D 1 » 

5070 

.?/ 

. 4737491 

74 

27320 

4599D1 , 

5030 

. 657925 1 

21 

201 01 

0099D1 ? 

5090 

.8525161 

*2'  *2' 

10654 

1 43 OD 1 , 

5 1  00 

o 

. 1060460 

\>  *7/ 

02745 

2502D2 , 

5 1  1 0 

y . 

. 1230132 

74 

so  OS  1 

46 9602? 

5120 

?/ 

. 1510441 

•7  cr 

73075 

5 1 53D2 , 

5  3  30 

?■/ 

. 1 750230 

7  0 

4  5  S' 7  3 

885802, 

5140 

Zi 

. 1993721 

44 

9566  1 

S362D2 , 

5  1  50 

9'. 

. 22552 1 6 

■2'  O 

53123 

4229D2 , 

5 1 60 

9. 

. 2519122 

1  1 

i~i 

68 1 5D2 > 

5170 

9. 

. 2739927 

1 3 

S3 840 

3916D2, 

5 1 30 

<>'. 

. 3067 1 36 

01 

06080 

672SD2 , 

5 1  90 

9, 

. 3350507 

34 

50 1 36 

83S9D2 , 

5200 

& 

. 3639544 

52 

0S033 

053602, 

5210 

•  6  V  -2'  o'  V  o  o 

41 

87 1 99 

4940D2 , 

5220 

e  _ 

. 423356 1 

6460753 

485002, 

5230 

?/ 

. 45330 1 3 

O 

93476 

90S0D2  , 

5240 

. 4347113 

13 

51335 

223902 , 

5250 

. 5 1 60667 

55 

67764 

3  736' D  2  , 

5260 

9, 

. 5473472 

■TV  O 

98 1  1 2 

3 1 92D2 . 

5270 

. 5300360 

52 

22980 

5 1 9902 , 

5230 

& 

. 6 1261 70 

17 

6 1 002 

002002 , 

5290 

9. 

. 6455753 

L 

»_■ 

27006 

33 1 1 02 , 

s-**** 
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rn  co 


3  1 0 

S<  .712570 3 S 9 6  7 1 630 0 9 0 D 2 

1 

320 

St  .  7 4353236343330 1 S43D2 

330 

/ 

340 

DPQ  =  DBLE (IP) 

350 

IFdP.LE.30)  STIRL  = 

LNFACT ( IP) 

I F ( I P . GT . 30 )  STIRL  = 

(  DPQ  +•  .  5  )  *  DLOG  (DPQ)  - 

370 

St  DPQ  +  K1  +  K2/DPC  -  K 

2/ ( . 3D2  *  DPQ  *  DPQ  *  DPQ) 

330 

RETURN 

390 

END 

340 OC - 

34  1 0 
5420 
5430 
5440 
5450 

54-30  4040 

5470 

5430 

5490 

5500 

55 IOC - 

3520 
33  30 
5540 
5330 
5530 

5530 

5530 

5600 

3-6  1  OC 

562'OC  **•*•*• 

5630C 

5640 

5650C 

3  6  0  C 

5670C 

563 OC 

5690 

5700C 


SUBROUT I NE  ■  PC I SSY ( P , R , E , FO , I STAR , IF I N ) 
DOUBLE  PREC I S I ON  P ,  R  >  E, P 0 
DO  4040  II  =  I ST AR+ 1,1 FIN 
R  =  R  -  DEXP(P) 

E  =  E  -  R 

P«P-  DLOG(DBLE< I  I ) >  +  PO 
IF  (E  .  LE .  0.000001)  E  0.0 
IF  (  R  .LE.  0.0  )  E  =  0.0 

RETURN 

END 


SUBROUTINE  S I ZOP ( I ) 

COMMON/ AVA/AVLOG ( 70 ) , CRANK ( 70 ) , IPS ( 70 ) 

COMMON./ BASE/BEBO 

COMMON/ NEXT/B80L I B  ? DBOLIB. PROLIB, DRBLIB 
COMMON/PARAMS/ I DONE , MAXRED , MODE • NCOMP , NOB , SPB 
COMMON/PR  I C £/ I COST ( 70 ) , NBEST ( 70 ) 

REAL  BEBO ( 7 0 ) 

DOUBLE  PREC  I  S I  ON  BBOL  IB ( IS >  »  DBOL  IB( 13) ,  F  ROL I B  (  1 S  >  ,  DRBL 1 3  <  1 S.) 
REAL  CRANK (IS) 

DETERMINE  THE  RANKING  VALUE  OF  EACH  PERMISSIBLE  PURCHASE 
DO  4545  N  =  NOB  -  MAXRED,  NOB 

*****THE  RANKING  VALUE  IS  THE  RATIO  OF  EBO  REDUCTION  TO  COST 
IF  AVAILABILITY  IS  UNDEFINED. 

IF  (  MODE  . EQ.  1  )  CRANK( N)  =  (BEBO(I)  -  BBOL I B ( N) ) /DBLE ( N ) 


57  IOC 
572 OC 
5730C 
5740  4545 
5750 


£t. 


****THE  RANKING  VALUE 
-  INCREASE  TO  COST, 


IS  THE  RATIO  OF  THE  LOG  OF  AVAILABILITY 
IF  AVAILABILITY  IS  DEFINED. 


IF  (  MODE  . EQ.  2  ) 

CRANK (N)  =  (  DBLE ( IPS( I ) )  *  DLOG (  1 . ODO  -  BBOL IB ( N ) / 

<  DBLE  <  I  PS  (  I  )  )  *DBLE  ( SP  B  )  ) .)  -  AVLOG  (  I)  )/ 

( DBLE ( N )  *  DBLE ( I  COST ( I ) )  ) 

IF  (  MODE  .EC.  1  )  CRANK (NOB+1)  =  (  BEBO ( I )  -  BBOLIB ( NOB+1 )  ) 

IF  (  MODE  .EC.  2  ) 

CRANK (NOB+1)  =  (  DBLE ( I PS ( I ) )  *  DLOG (  1 . ODO  -  BBOLIB < NOB+ 1 > / 

( DBLE ( I PS ( I ) ) *DBLE ( SPB ) ) )  —  AVLOG ( I ) ) / 

DBLE ( I COST ( I ) ) 
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5840*:: 

3850C 
5860C 
5870 
5880 
5890 
5900 
59 1 0 
5920 
3920 
5940 
5950 

5960C - - - 

5970  SUBROUTINE  COMPOP ( I MAX ) 

5980  COMMON/ AVA/AVLOG ( 70 ) , CRANK ( 70 ) , I PS ( 70 ) 

5990  COMMON / P ARAMS / 1 DONE , MAXRED? MODE , NCOMP , NOB,  SPB 

6000  REAL  TEMP (70) 

60 1 0  I NTEGER  I POS f  70 ) 

6020 C 

THIS  IS  THE  FIRST  RANKING  OF  THE  RUN,  RANK  ALL  OF 
COMPONENTS . 


*****NOW  CHOOSE  THE  MOST  COST  EFFECTIVE  PURCHASE, 

CRANK ( I )  -  QRANK ( NOB+ 1 ) 

NBEST(I)  =  NOB  +  1 

DO  5545  N  =  NOB  -  MAXRED,  NOB  +  1 

I F (  QRANK ( N )  . LE -  CRANK ( I )  )  GO  TO  5545 

CRANK ( I )  =  QRANK (N) 

NBEST(I)  =  N 
5545  CONTINUE 
RETURN 
END 


IF  (  I DONE  . EQ-  1  >  GO  TO  5550 

***** TEMPORARILY  STORE  THE  COMPONENT  RANKING  VECTOR  3 

DO  5150  L  =  1 , NCOMP 
TEMP ( L )  =  CRANK ( L  > 


6U0UL 
60401 
6050*: 

6060 

60  TOC 
60S0C 
6090C 

6 1  00 

6110  5150 
6 1 20C 

6130C*****RANK  EACH  COMPONENT.  ONCE  IT  IS  RANKED,  SET  ITS  TEMP* 
6 1 40 C  RANKING  VALUE  TO  ZERO. 

6 1  50  C 
6 1 60 
6  1 700 
6130C 
6 1 90C 
6200 

62  IOC 

6220C 
6 2 4 0 C 
625  OC 
6260C 
6  2 700 
62S0  5250 
6290C 

6300  5S50  TEMP ( IPOS ( M)  )  =  0 -  0 
63 1  OC 

6320  GO  TO  5850 

6330C*****UPDATE  RANKINGS  TO  REFLECT  THE  LATEST  PURCHASE 
6340C 
6  3  50 1 

6360C»****FIND  THE  NEW  RANKING  POSITION  OF  THE  MOST  RECENTLY  PUI 
ART  or 


HE 


T  E  M  F' . 


RAF:' 


DO  5350  M  =  1 , NCOMP 

*****HYPOTHESISs  COMPONENT  1  SHOULD  BE  RANKED  IN  POS  I  I  -t  ON  M. 
I  POS:  <  M  )  =  1 
DO  5250  L  =  1  •>  NCOMP 


*****TEST  HYPOTHESIS:  IS  ANY  UNRANKED  COMPONENT  MORE 

EFFECTIVE  THAN  THE  ONE  CURRENTLY  RANKED  IN  POSIT 
IF  SO,  ASSIGN  THE  MORE  EFFECTIVE  COMPONENT  TO  M. 

IF  (  TEMP ( L )  . OT.  TEMP ( I POS ( M ) )  )  IPOS(M)  =  L 


COST 
ON  M? 


CHASED 
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6.330  3330 

NCWP.OS  =  1 

639'J  5650 

NN  =  NEWPOS  -i-  1 

6400C 

64  1  OC 

*****SEE  IF  THIS  VALUE  OF  NEWPOS 

6420C 

IF  IMAX. 

6430C 

6440 

IF  (  CRANK (I MAX)  .GE.  CRANK 

6450C 

6460C 

***-**IF  IT  IS  NOT,  THEN  THE  COMP 

647 OC 

BUMPED  UP  ONE  POSITION. 

6460C 

6490 

I F'OS  <  NEWPOS  )  =  I POS  (  NN  ) 

6500 

NEWPOS  =  NN 

63  1 0 

IF  <  NN  . LT.  NCOMP  )  GO  TO 

6520C 

6330C 

**■*•«••* IF  IT  IS,  THEN  RANK  IT. 

6540C 

6550  5750 

I FOS( NEWPOS)  =  IMAX 

■656 OC 

6570  5 S' 50 

IMAX  =  I POSC 1) 

65S0 

RETURN 

6590 

END 

66  1 0 

SUBROUTINE  UPDATE ( I ) 

6 6 2 0 C  C OMM 0 N /  E B 0 N /  BEBON  (  7 0  )  ,  DE B 0 N  (70) 

6630  C0MM0N/PARAM3/ID0NE, MAXRED, MODE , NCOMP ,  MOB ,  SPB 

6640  COMMON/PR I CE/ I  COST ( 70 ) , NBEST < 70 ) 

6630  COMMON / SPARE / *JSP ARE  ( 70  )  ,  K SPARE  (  70  ) 

6660  COMMON / S  YST / S YSLOO , SYSTAV , TOT C ST 

66 700 

66S0C*****UPDATE  THE  SPARES  INVENTORIES  AND  THE  TOTAL  COST  COUNTER. 
6690C 

3EB0  < I )  =  BEBON ( I ) 

DEBO(I)  =  BEBON < I ) 

■JSPARE  (  I  >  =  J SPARE  (  I  >  -  NOB  +  NBEST  (  I  ) 

I F (  NBEST ( I )  .  EQ.  NOB  +  1  )  GO  TO  4433 

KSF'ARE (  I  )  =  KSPARE(I)  +  1 
TO TOST  =  TOTCST  +  NBEST ( I )  *  ICOST(I) 

4433  TOTCST  =  TOTCST  +  ICOST(I) 

67S0C 

6790C*****UPDATE  THE  PURCHASED  COMPONENT'S  AVAILABILITY.  ALSO  UPDATE 
6S00C  THE  SYSTEM'S  AVAILABILITY  IF  IT  IS  COMPUTABLE. 

63 1 00 

6320  3333  CALL  AVAIL(I) 

6330  RETURN 

6340  END 

63500 - 

6360 

6330 
6390 
6900 
69  i  o 


6700C 

67 1 OC 

6720 

6730 

6740 

6730 

6.760 

6770 


SUBROUTINE  AVAIL(I) 

COMMON/ AVA/AVLOG ( 70 ) , CRANK ( 70 ) , IPS ( 70 ) 

:  0  M  M  0  N  /  B  A  S  E  /  B  EB  0 

COMMON/PARAMS/ I DONE , MAXRED , MODE , NCOMP , NOB , SPB 
COMMON/PRICE/ I  COST ( 70 ) , NBEST ( 70 ) 

COMMON / S YST / SYSLOO , SYSTAV, TOTCST 
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6920 
6930i 
694< 

6930C 

6  96  Of 
697' 
6930 
6990 


REAL  BE DO < 70 ) 


s-****C0MPUTE  THE  LN  OF  THE  COMPONENT  AVAILABILITY.  IF  THE  COMP 
HAS  ALREADY  BEEN  RANKED ?  USE  THE  QUICK  METHOD  OF  COMPUTA¬ 
TION,  STATEMENT  NUMBER  6660. 


IF  <  I  DONE 


EQ . 


MODE  . EQ.  2  )  GO  TO  6660 


IF  (  BEBG(I)  .  GE.  FLOAT  ( IF’S(  I  )  )  *  SF'B  )  GO  TO  9960 


7000 

AVLOG(I)  =  DBLE ( I PS ( I ) )  *  DLOG< 

1  .  ODO  —  BEBO  f  .1  )  / 

7010  S< 

( DBLE ( I  PS ( I ) ) *DBLE ( SP 

B  >  )  ) 

7020 

SYSLOC  =  SYSLOG  +  AVLOG ( I ) 

7030 

MODE  =  2 

7040 

GO  TO  9060 

7050  6660 

IF  (  NBEST(I)  .LE.  NOB  ) 

7060  & 

DELTA  =  DBLE ( NBEST < I ) >  *  DBLE (I COST 

(  I  ) )  *  OBLE ( OR AN I 

7070 

IF  (  NBEST ( I )  .EQ.  NOB  +.1  ) 

7030  S< 

DELTA  =  DBLE '(I  COST  (  I)  )  *  DBLE  <  CRANK 

(  I  )  ) 

7090 

AVLOG ( I )  =  AVLOG  < I )  +  DELTA 

7 1  00 

SYSLOG  =  SYSLOG  +  DELTA 

7 1  1 0 

IF  (  SYSLOG  . GE.  -Sl.ODO  )  SYSTAV  = 

EXP ( SYSLOG ) 

7120  99 60 

CONTINUE 

7 1 30 

RETURN 

7 1  40 

END 

7 1 50  C 
7 1 60 
7170 
7 1  SO 
7  1  90 
7200 
7210 

7230 

7240 

7250 

7270 

72301] 

7290i] 

7300C 

73  ICC 

7 320 C 

733 OC 

7 340 L 

7350 

7360 

7370 

733 OC 

7390C 

7400C 

7410 

7^20 

7430 

7440 

74=;o 


JBROUTINE  PRELIM  a) 

OMMON/ EASE/ BEBO 

OMMON /DEPOT  /DEBO,  DMLOG?  DP  ROB?  DREBO 
OMMON/ MEAN /BMEANO ( 70 ) , DMEAN ( 70 ) 

OMMON /OPTS/ IPRE 

OMMON /P ARAMS/ I DONE , MAX RED , MODE ? NCOMP , NOB ? 3PB 


OMMON / PR I CE / I C 
OMMON / SPARE / US 
OMMON/SYST /SYS 


;OST  (  70  )  ,  NBEST  (  70  ) 

JPARE  (  70 )  ,  K SPARE  (  70  ) 

JLOG?  SYSTAV,  T0TC3T 
DOUBLE  PR EC 1 3 1 ON  AM , PO ? P , R ? U 
REAL  BEBO ( 70 ) 

REAL  DEBO ( 70 ) , DM LOO ( 70 ) , DP ROB ( 70  > , DREBO ( 70 ) 


******  MAKE  PRELIMINARY  BUYS  OF  EACH  COMPONENT  FOR  THE  DEPOT 
THE  BASES,  IF  NECESSARY. 


*****TEST:  DOES  THE  DEPOT  NEED  A  BUY? 

IF  <  (  IPRE  .  EQ.  1  .AND.  DMEAN {  I  >  LE.  10.0)  .OR. 
Sc  (  I  F’RE  .  EQ .  2  .  AND .  DMEAN  (  I  )  .  LE .  33 . 0 )  ) 

GO  TO  1005 

******  IF  IT  DOES  i  SET  THE  DEBO  VARIABLES  AND  BUY .  ***** 

P  =  DBLE  ( DF'ROBC  I  )  ) 

R  =  DBLE ( DREBO < I ) > 

E  =  DBLE ( DEBO < I ) ) 

PO  =  DBLE  (  DMLOG  (  I  )  ) 

AM  =  DBLE ( DMEAN < I ) ) 
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7  4  AO 

7470C 

743 OC 

7490C 

7300 

75 1  OC 

7520C 

7330C 

7540C 

7330 

756  OC 

737 OC 

7530 C 

7590 

7600 

7610 

7620 

7630C 

7640C 

7630C 

7660 


7330C 
7340 C 
7350 C 
736 OC 
7870 
7  8  S'OC 
7890C 
7900 C 
79 1 0 
7920 
7930 
7940 
7950 


LLL  =  IF1X(DMEAN(I)> 

*■**•*■* ;[F  USER  OPTED  FOR  2,  BUY  A  PIPELINE  S  WORTH. 

I F  (  IPPE  .  EQ .  2  )  CALL  LI B  <  P ,  R ,  E ,  PO ,  AM ,  .JSF'ARE  (  I  )  ) 

«■**** IF  USER  OPTED  FOR  1,  BUY  SIX  STANDARD  DEVIATIONS 
LESS  THAN  A  PIPELINE'S  WORTH. 

IF  (  I PRE  .  EQ .  1  )  CALL  CON  (  P ,  E ,  PO ,  AM ,  JSF'ARE  (  I  )  ,  LLL  ) 

ACCOUNT  FOR  THE  GLUMF'  AT  THE  DEPOT 

DP ROB ( I )  =  P 
DREBO ( I )  =  R 
DEBO ( I )  =  E 

TOTCST  =  TOTCST  +  JSPARE(I)  *  lCOST(I) 
*****COMPUTE  BMEAN  GIVEN  THE  PURCHASED  SPARES.**** 
AM  =  DBLE < B ME A NO < I ) )  -  < AM-E ) /DBLE ( NOB ) 


7670 

E 

=  AM 

7  6  S'OC 

7690C 
770 OC 

*-*■*-**  TE 

ST :  DOES  THE 

BA: 

BE  NEED  A  BUY? 

77  1 0 

IF 

(  (  1 PRE  . EQ 

«  4 

.AND.  AM  . LE.  10.0  ) 

.  OR . 

7720 

& 

(  I PRE  .EQ 

.  1 

.  AND .  AM  .  Li_ .  .  O  ? 

) 

7730 
7740 C 

GO 

TO  1015 

7760C 

-X.  -X.  DC  -X.  -X.  T  p" 

IT  DOES 7  BE 

T  T 

HE  BEBO  VARIABLES  AND 

w*.:|  jy  _X.  DC  DC  DC  _X. 

777 0 

PO 

=  DLOG(AM) 

7790C 

•  LLL 

=  I FIX (AM) 

7300C 

■K- I F 

THE  USER  OPTED 

FOR  2,  BUY  A  PIPLINE' 

S  WORTH. 

73 1  OC 
7320 

IF 

(  I PRE  .EQ. 

2  ) 

CALL  LIB<P»R»E» PO ? AM . 

KSPARE (  I  ) 

IF  THE  USER  OPTED  FOR  1,  BUY  SIX  STANDARD  DEVIATIONS 
than  A  PIPELINE'S  WORTH. 

I F  (  I PRE  . EQ .  1  )  CALL  CON ( F, £ , PO » AM , KSPARE <  I  )  , LLL ) 

*****UFDATE  BEBO  AND  COST***** 

1015 


BEBO (I)  =  E 
TOTCST  =  TOTCST 


KSPARE ( I )  *  I COST (I)  *  NOS 


1005  CONTINUE 

RETURN 
END 


TIME  OUT 
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EPILOG 


1  1st 


I  COOL 
1 0 1 0 
1 020 
1 030 
1 040 
1 050 
1060 
1 070 

1 08  OC 
1 090C 

I I  ooc 


JTF'UT .  ****■! 


*****THIS  PROGRAM  PROCESSES  DRAMA  OU' 

I NT EGER  K SPA RE ( 70 ) , J SPARE < 70 ) 

I NTEGER  I PS ( 70 ) , I COST ( 70 ) 

IF  (  NAVA I L  . EQ .  0  )  NT ABLE  =  27 
IF  (  NAVA I L  . GT.  0  )  NT ABLE  =  NAVA I L 
REAL  EEBO ( 70 ) , DEBO ( 70 ) 

CALL  FMED I A  <  9 »  0 ) 

CALL  FMEDI A ( 10,0) 

ALL  FMED IA( 11,0) 

ALL  FMEDI A (12,0) 

ALL  FMEDI A (13,0) 


1  1  IOC 

U20C*****USER  OPTION  VARIABLES.  ***** 
1130  I  OUT  =  1 

1140  NAVA I L  =  0 

1 1 50  C 


1 1 60C 

INPUT  SCEN 

1 170 

READ (9, 10 

1 1 30C 

***** input 

1 1 90 

IF  (  I  OUT 

1200 

DO  999  I 

1210 

■“/■yy 

READ (9,  10 

1  LLC) 

1  09 

FORMAT ( V ) 

1230 

1 000 

CONTINUE 

1240C 
1 25  OC 
1 2600 
1270C 
1230 
1290 
1 300 
1310 
1320 
1 330 
1  34  OC 
1 350C 


*****PRINT  AN  OUTPUT 
USER 


TABLE  FOR  EACH  AVAILABILITY  SELECTED  BY  TH 


DO  8000  IT  =  1 , NT ABLE 
READ (10)  K SPARE 

IF  (  I  OUT  . GT.  1  )  READ (11) JSPARE 

IF  (  I  OUT  . GT.  2  )  READ (12)  BEBO 

Ic  (  I OUT  . GT.  2  )  READ (13)  DEBO 


IF  (  NAVA  I L  .LT.  NT ABLE 


AND.  NAVA  I L  . GT.  0  )  GO  TO  3000 


***** INITIALIZE  COUNTERS  IF  I  OUT  = 


-X.  JJ,  .X.  -X. 


1 360 

IF  (  I  OUT 

.  LT. 

3 

) 

GO  TO  2000 

1370 

TCOS- T  =  0. 

0 

1 330 

DCOST  =  0 . 

o 

1 390 

BOOST  =  0 . 

o 

1 400 

DBG TOT  =  0 

.  0 

1410 

BBOTOT  =  0 

.  0 

1420 

AVTOT  =  1 . 

0 

1430  2000 

CO NT I NUE 

1 440C 

1450C 

*****PRINT 

A  TABLE 

HEADER. ***** 

1 460 

IF  (  I OUT 

.  EQ. 

1 

) 

WRITE  (6,3016)  IT 

1470 

IF  (  I  OUT 

.  EQ. 

2- 

) 

WRITE (6, 3026)  IT 

1 430 

IF  (  I OUT 

.  GT. 

L 

> 

WRITE (6, 3036)  IT 

1490  8016 

FORMAT ( 1 H 1 

,  20  (  •" 

-X- 

) 

, -BASE  STOCKAGE  LEVELS 

.1 500 

?<  ■-  LEVEL  •-  , 

1 3 ,  1  1  (  •" 

•-  )  ,  /  /  ) 

1510  3026 

FORMAT ( 1 H 1 

,  i  o  ( 

-X. 

) 

,  •■’TOTAL  STOCKAGE  LEVEL 

II-  T  I  T  \  / 
.  M"!  1  _  1  *  7 
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*****OUTPUT  DATA  TO  TERMINAL.  PERFORM  NECESSARY  COMPUTATIC 
IF  OPTION  3  OR  4  HAS  BEEN  SELECTED. 

DO  7000  1=1, NCOMP 
I TAB  =  0 

#****3TART  A  NEW  PAGE  EVERY  40  L I NES .  ***** 

IF (  MOD (10,1)  .EQ.  1  )  ITAB  =  1 

IF  (  ITAB  .EQ.  1  .AND.  I  . GT.  1  )  WR I TE ( 6 , 7006 ) 

7006  FORMAT <1 HI) 

IF  (  I  OUT  .EQ.  2  )  GO  TO  6200 
IF  (  I OUT  . GT.  2  )  GO  TO  6300 
*****OUTPUT  FOR  OPTION  1 . ***** 

IF  (  ITAB  . EQ.  1  )  WRITE (6, 7016) 

70 1 6  FORMAT  ('  "‘  ,  •"  COMPONENT  TYPE  "  ,  SX  ,  BASE  STOCK •" ,  / /  ) 

WR I TE ( 6 , 70 17)  I ,  KSP ARE ( I ) 

7017  FORMAT  (  •"  •" ,  1  1  X  ,  13,  1 1  X ,  1 4  ) 

GO  TO  7000 

***-**OUTPUT  FOR  OPTION  TWO .  ***** 

6200  IF  (  ITAB  .EQ.  1)  WR I TE ( 6 , 7026 ) 

WR I TE  <  6 ,  7027  )  I , USFARE ( I ) , KSP ARE ( I ) 

7026  FORMAT (■-  •"  ,  •" COMPONENT  TYPE  ,  5X,  'DEPOT  STOCK  " ,  3X,  '  BASE 
Se  •"  STOCK",//) 

7027  FORMAT  (  "  •" ,  1 1  X  ,  13,  1 2X  ,  1 4,  1 1  X  ,  1 4 ) 

GO  TO  7000 


:s  FOR", 


DNS 


1 S 1  oc 

1320  *****OUTF'UT  FOR  OPTIONS  3  AND  4***** 

1330  6300  CONTINUE 


1 3400: 

*  *  *  *  <•  C  OM  F  U  T  E 

COSTS  AND  AVAILABILITIES***** 

1350 

LCD? 

T  =  (  NOB*KS 

PARE  (  I  )  +  JSF'ARE  (  I  )  )  *  ICOST(I) 

1360 

T  COS 

T  =  TOO ST  + 

LCOST 

1370 

DCOS 

T  =  DCOST  + 

JSF'  ARE  (  I  )  *  I  COST  (  I  ) 

1 330 

BCOS 

T  =  BOOST  + 

( NOB*K SPARE ( I )*I COST  < 3 ) ) 

1 390 

AVCC 

MR  =  (  1 . 0  - 

BEBO  (  I  )  /  (  FLOAT  (  I  PS  (  I  )  )  *  SPB 

1 900 

it  **  I  PS  (  I  ) 

1 9 1 0 

AVTOT  =  AVTOT  * 

AVCOMP 

1  920 

BBOTOT  =  BBOTOT 

+  BEBO  ('  I  ) 

1 930 

DBOTOT  =  D 80 TOT 

+  DEBO ( I ) 

1 940C 

*  -x-  *  *  *  0  UTF'UT 

COMPONENT  DATA***** 

1  930 

IF  ( 

ITAB  . EQ.  1 

)  WRITE (6, 7036) 

1  960 
1 970 
1 930 
1990 
2000 
20 1 0 
2020 
2030C 
2040C 
2050 


WRITE (  6,7037  ) 

?■:  I ,  .J  SPARE  (  I  )  ,  KSP  APE  (  I  )  ,  LCOST,  DEBO  (  I  )  ,  BEBO  (  I  )  ,  AVCOMP 

FORMAT  (  •"  •"  ,  •"  COMP .  TYPE "' ,  3 X' ,  "  DEPOT  STOCK  •" ,  3 X  ,  •"  BASE  STOCK  "' , 
?•:  3 X  ,  "'  BASE  EBOi ,  3 X  ,  •"  DEPOT  EBO  -  ,  3 X  ,  •"  BASE  EBO "' ,  3 X  , 

S<  "'  A VA 1  LAB  I L I T Y "',//) 

7037  FORMAT  (  •"  •"  ,  7X  ,  13,  10X  ,14, 9X  ,  13,  3X  ,  F9. 3, 2X  ,  F9. 3,  5X  ,  F10. 6  ) 
7000  CONTINUE 


7( 


IF  (  I  OUT  . LT. 


)  GO  TO  SOOO 
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20600 

2070C*****WRITE  A  S 

20SOC 

2090 

WR I TE  < 6 , 

2100  9006 

FORMAT ( ' 

2110  2< 

13,///) 

2120 

WRITE (6, 

2 1 30 

WR I TE ( 6 , 

2140 

WR I TE ( 6 , 

2 1 50 

WRITE (6, 

2160 

WRITE (6, 

2170 

WRITE (6, 

2 ISO  90 1 6 

FORMAT  < ' 

2190  9026 

FORMAT ( ' 

2200  9036 

FORMAT ( ' 

22 1 0  9046 

FORMAT ( ' 

2220  9056 

FORMAT  (  •- 

2230  9066 

FORMAT  (  •- 

2240  3000 

CONTINUE 

V  !“i{_) 

STOP 

2260 

END 

■K- 


UMMARY  OUTPUT  TABLE  IF  ONE  13  CALLED  FOR.  #*#■*# 
9006)  IT 

*  1 H 1 7  •"  END  I TEM  SUMMARY  FOR  AVA I  LAB  I L I TY  LEVEL  ' , 

9016)  BOOST 
9026 )  BOOST 
9036)  TCOST 
9046 )  DBOTOT 
9 0 56)  B  B  0 T 0 T 
9066)  AVTOT 

"A  "TOTAL  COST  OF  DEPOT  SPARES' .24X, FI 0. 0) 

•-  ,  -  TOTAL  COST  OF  BASE  SPARES  '  ,  25X  ,  F 1 0 . 0  ) 

,  -  GRAND  TuTAL  OUST  , 34X , F10. 0 ) 

'.'REMAINING  EBO  AT  DEPOT', 24 X , FS.  3 ) 

'  •' ,  •-  REM  A I N I  NO  EBO  AT  EACH  BASE  '  ,  20  X  ,  FS .  3  ) 

' ,  ' END- 1 TEM  AVA I  LAB I L I TY ' , 26 X , F 1 0 . 6 ) 
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APPENDIX  C 


EXECUTING  JOBS  ON  THE  HONEYWELL  635 


TIME -SHARING  MODE 

When  DRAMA  is  being  run  for  a  small  data  base  (250  items  or  less)  with 
moderate  demand  (an  average  mean  base  demand  across  all  components  of  10  or 
less)  and  for  a  small  number  of  bases  (15  or  less),  it  is  practicable  to  run 
the  model  in  a  time-sharing  mode.  Setting  up  a  time-sharing  run  is  quite 
simple,  for  all  that  needs  to  be  done  (besides  setting  up  the  input  data  base 
with  PROLOG)  is  to  "attach"  the  input  and  output  files  to  the  main  program. 
The  easiest  way  to  do  this  is  to  add  a  line  of  code  to  the  main  program  for 
each  file  that  needs  to  be  attached.  For  example,  the  code  required  to  attach 
the  input  file  XM-1DATA  is 

CALL  ATTACH(9,MOS29/N232D/XM-lDATA;n,  3,0,ISTAT,) 

The  9  in  the  parentheses  is  the  logical  file  number  for  DRAMA  inputs. 
The  string  OS29  is  a  catalog  string  (which  differs  among  using  organizations) . 
The  string  N232D  is  a  project-number  string  (which  differs  among  users  within 
a  using  organization) . 

If  one  wishes  to  see  additional  outputs  he  must  attach  the  appropriate 
output  files.  The  appropriate  logical  file  number  for  each  type  of  output 
file  is  designated  below: 

File  Type  Logical  File  Number 

Base  Spares  10 

Depot  Spares  11 

BEBO  12 

DEBO  13 
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Therefore,  if  one  wishes  to  output  the  Base  Spares  data  to  a  file  called 
BSPARES,  he  must  add  the  following  statement: 

CALL  ATTACH (9 , "0S29/N232D/BSPARES ; " ,  3,0,  ISTAT,) 

RUNNING  IN  REMOTE  BATCH  MODE 

Jobs  too  large  for  time-sharing  can  usually  be  run  in  remote  batch  mode. 
To  run  in  remote  batch  the  user  must:  (l)  put  the  DRAMA  program  into  a  form 
submittable  in  remote  batch,  (2)  properly  define  the  input  and  output  files  in 
a  JCL  program  that  will  submit  DRAMA  in  remote  batch,  and  (3)  establish 
suitable  limits  in  the  JCL  program  for  DRAMA'S  running  time,  core  requirement, 
and  output  requirement. 

To  put  DRAMA  into  a  form  that  can  be  submitted  in  remote  batch,  the  user 
must  change  all  CALL  ATTACH  statements  in  the  program  to  CALL  FMEDIA 
statements.  For  example,  he  must  change  the  statement  that  attaches  the  input 
file,  XM-1DATA,  from  , 

CALL  ATTACH  (9 , "0S29/N232D/XM-1DATA; " ,3 , 0 , ISTAT, ) 
to 

CALL  FMEDIA  (9,0) . 

Similarly,  he  must  change  the  statement  that  attaches  the  output  file 
BSPARES  from 

CALL  ATTACH  (10 , "0S29/N232D/BSPARES ; " ,3 , 0 , ISTAT, ) 
to 

CALL  FMEDIA  (10,0). 

To  make  sure  that  the  correct  file  is  input  as  logical  file  9  (and  to 
make  sure  that  the  correct  files  are  output  as  logical  files  10,  11,  12  and 
13),  the  user  must  include  a  PERMFILE  or  PRMFL  statement  in  his  JCL.  For 
example,  to  input  the  file  XM-1DATA,  the  user  should  employ  the  statement 

$ : PRMFL : 09 , R , S , 0S29/N232D/XM- 1DATA 
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In  this  statement,  09  is  the  logical  file  number  that  will  represent  the 
file  XM-1DATA  in  the  main  program,  the  letter  R  indicates  that  this  file  will 
be  read,  the  letter  S  indicates  that  it  is  sequential,  and  the  string 
0S29/N232D/XM- 1DATA  is  the  complete  catalog/project  number/file  name  of  the 
input  file. 

To  assign  the  output  file  MBSPARESff  to  the  correct  logical  file  number, 
the  following  is  required: 

$ :  PRMFL: 10 ,W,S ,0S79/N232D/BSPARES 

The  number  10  here  stands  for  the  logical  file  number  of  the  output  file, 
and  the  letter  W  indicates  that  this  logical  file  must  be  written  to,  not 
read.  For  each  logical  file  in  the  main  program,  there  must  be  a  PERMFILE 
statement  in  the  Job  Control  Language  program. 

LIMITS  FOR  REMOTE  BATCH  JOBS 

The  limits  for  a  remote  batch  job  will  need  to  be  altered  with  the  size 

of  the  job.  The  running-time  limit  of  the  job  is  driven  by  the  number  of 

components  and  by  the  size  of  BMEAN0  and  DMEAN.  For  a  data  base  of  200 
components  in  which  the  average  value  of  BMEAN0  was  less  than  10,  20  CPU 
seconds  were  found  to  be  an  adequate  limit.  When  average  demand  is  higher-- 
say  25 --then  a  100-CPU-seconds  limit  would  be  more  reasonable. 

Core  requirements  are  driven  primarily  by  the  number  of  components.  Data 

bases  of  200  components  or  less  should  require  no  more  than  15K  of  core. 

Larger  data  bases--say  1000  components--will  probably  require  more  than  30K  of 
core. 

Output  requirements  of  DRAMA  are  minimal,  so  we  have  arbitrarily  estab¬ 
lished  a  nominal  number  of  lines--2K--as  the  DRAMA  output  limit. 

If  a  job  exceeds  one  of  its  limits  it  will  terminate,  and  the  user  will 
have  to  resubmit  it  with  altered  limits  in  order  to  obtain  a  complete  run. 
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The  LIMITS  statement  of  the  Honeywell  JCL  looks  like  this: 

LIMITS: 20, 32K, ,2K 

In  this  statement  20  is  the  time  limit  in  CPU  seconds,  32K  is  the  core 
limit  in  6-bit  words,  and  2K  is  the  output  limit  in  lines. 

THE  JCL  PROGRAM 

A  sample  JCL  program  for  a  complete  DRAMA  run  is  presented  below: 

10##N, J 

20$ : IDENT :OS2033N232D  ,OS29UFABBRO 

30$:LIMITS:20,20K, , 4K 

6  0  $ : OPT I ON : FORTRAN 

70$: FORTY :NDECK 

90$ : SELECTA:OS29/N232D/DRAMA 

95$ : LIMITS : 2 , 20K, , 2K 

100$: EXECUTE 

110$:LIMITS:20,20K, , 2K 

120$ : PRMFL : 09 , R, S , OS29/N232D/XM-1DATA 

125$ : PRMFL : 10 , W, S , OS29/N232D/BSPARES 

140$: END JOB 

The  IDENT  statement  contains  the  user's  complete  catalog-string.  The 
first  LIMITS  statement  defines  the  limits  of  the  entire  job;  the  OPTION  state¬ 
ment  denotes  the  language  of  the  submitted  program;  the  FORTY  statement 
instructs  the  FORTRAN  compiler  to  read  the  stored  program  as  if  it  were  a  deck 
of  cards;  the  SELECTA  statement  gets  DRAMA;  the  second  LIMITS  statement  estab¬ 
lishes  the  limits  for  the  compilation  of  DRAMA;  the  last  LIMITS  statement 
establishes  limits  for  DRAMA'S  execution;  and  the  PRMFL  statements  define  the 
logical  input  and  output  files. 

RUNNING  THE  JCL  PROGRAM 

To  run  the  JCL  program,  the  user  must  enter  the  Honeywell's  CARDIN  sub¬ 
system.  He  does  this  by  keying  in  "CARDIN"  when  the  operating  system  asks  him 
what  subsystem  he  wishes  to  use. 

Once  in  CARDIN,  the  user  can  run  the  program  by  putting  the  JCL  program 
into  his  workspace  and  then  typing  RUN.  The  operating  system  will  respond  to 
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this  command  by  typing  out  a  job  reference  number,  or  SNUMB.  This  SNUMB  must 
later  be  used  to  retrieve  output. 

The  recommended  way  to  retrieve  output  is  to  type  in  the  command 
JOUT  "SNUMB". 

after  giving  the  program  a  reasonable  time  to  execute.  Then  the  user  can 
inspect  his  output  from  a  remote  terminal  by  using  the  Honeywell  Time-Sharing 
System/Batch  Interface  facility.  Finally,  to  get  hard  copy  of  his  output,  he 
should  type  in  the  command 
DIRECT  ONLINE 

which  will  cause  the  output  to  be  printed  out  at  the  line  printer  in  the  Air 
Force  Data  Services  Center.  Output  will  usually  be  printed  out  at  the  center 
within  a  few  hours  of  the  DIRECT  ONLINE  command. 
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